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hr  abstract  This  investigation  treats  the  problem  of  the  scattering  of  a  Rayleigh  wave 
bv  the  edge  of  a  thin  layer  which  covers  half  the  surface  of  an  elastic  half-space. 

The  ££ /action between  the  layer  and  the  half-space  is  described  approximatel^by 
means  of  a  model  in  which  the  effect  of  the  layer  is  represented  by  a  pair  of  boundary 
conditions  at  the  surface  of  the  half-space.  Two  parameters  -  one 

and  the  other,  stiffness  -  are  found  to  characterize  the  layer.  The  incident  Rayle  g 
wave  impinges  normally  upon  the  plated  region  from  the  unplated  side. 

In  the  case  where  the  mass  of  the  layer  vanishes,  the  problem  is  solved 
exact!  v  using  Fourier  transforms  and  the  Wiener-Hopf  technique,  and  numerical 
results  are  obtained  for  the  amplitudes  of  the  reflected  and  transmitted  surface .waves, 
i  In  the  more  general  case  of  a  layer  possessing  ooth  mass  and  stiffness,  a  perturba 
i  tion  procedure  leads  to  a  sequence  of  problems,  each  of  which  may  be  solved  using 
Fourier  transforms.  The  zeroth-  and  first-order  problems  are  solved  and ^the  re- 
suiting  approximate  reflection  and  transmission  coefficients  are  evaluated  numen 
ecliy  for  various  ratios  of  layer  mass  to  stiffness. 
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Abstract 

This  investigation  treats  the  problem  of  the  scattering  of  a 
Rayleigh  wave  by  the  edge  of  a  thin  layer  which  covers  half  the  sur¬ 
face  of  an  elastic  half- space.  The  interaction  between  the  layer  and 
the  half-space  is  described  approximately  by  means  of  a  model  in 
which  the  effect  of  the  layer  is  represented  by  a  pair  of  boundary 
conditions  at  the  surface  of  the  half-space.  Two  parameters  —  one 
representing  mass  and  the  other,  stiffness  —  are  found  to  character¬ 
ize  the  layer.  The  incident  Rayleigh  wave  impinges  normally  upon 
the  plated  region  from  the  unplated  side. 

In  the  case  where  the  mass  of  the  layer  vanishes,  the  problem 
is  solved  exactly  using  Fourier  transforms  and  the  Wiener-Hopf 
technique,  and  numerical  results  are  obtained  for  the  amplitudes  of 
the  reflected  and  transmitted  surface  waves.  In  the  more  general 
case  of  a  layer  possessing  both  mass  and  stiffness,  a  perturbation 
procedure  loads  to  a  sequence  of  problems,  each  of  which  may  be 
solved  using  Fourier  transforms.  The  zeroth-  and  first-order 
problems  are  solved  and  the  resulting  approximate  reflection  and 
transmission  coefficients  are  evaluated  numerically  for  various 
ratios  of  layer  mass  to  stiffness. 

*  This  investigation  was  carried  out  with  the  financial  assistance  of 
Hughes  Aircraft  Company  through  the  Howard  Hughes  Doctoral 
Fellowship  Program;  it  was  also  supported  in  part  by  the  Office  of 
Naval  Research  under  Contract  N0014-67- A-0094-0020. 
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The  present  work  is  concerned  with  the  scattering  of  an 
elastic  surface  wave  by  the  edge  of  a  thin  surface  layer. 

The  principal  motivation  for  studying  the  relation  between 
elastic  surface  waves  and  thin  .surface  films  comes  from  the  field  of 
microwave  signal  processing.  Electromagnetic  waves  may  be  con- 
verted  to  elastic  waves  in  solids,  and  vice  versa,  by  means  of  piezo- 
electric  transducers.  Operations  such  as  delaying  and  filtering  of 
signals  may  be  performed  on  the  elastic  waves  by  suitably  designed 
devices  [1]  .  The  advantage  of  carrying  out  such  operations  in  the 
elastic  -  rather  than  the  electromagnetic  -  medium  comes  about 
primarily  from  the  fact  that  the  phase  velocities  of  elastic  waves  are 
typically  five  orders  of  magnitude  smaller  than  those  of  electro¬ 
magnetic  waves  [2]  .  Thus,  for  example,  the  length  of  delay  line 
required  for  a  given  delay  is  smaller,  by  this  ratio  of  phase  veloc¬ 
ities,  for  an  elastic  wave  than  for  an  electromagnetic  wave. 

Elastic  surface  waves  prove  to  be  more  advantageous  than 
body  waves  because  of  their  accessibility.  Thus,  for  example,  sur¬ 
face  wave  delay  lines  may  be  built  with  numerous  taps  along  the  way, 
providing  incrementally  varying  delays. 

The  surface  of  an  elastic  solid  acts  in  a  sense  as  a  waveguide, 
i.  e.  ,  energy  once  transferred  into  surface  waves  tends  to  remain  in 
fields  concentrated  near  the  surface.  For  purposes  of  microwave 


Reference  fl]  provides 
of  various  surface  wave 


a  general  survey  of  the  design  and  function 
signal  processing  devices. 
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signal  processing  it  is  also  desirable  to  confine  —  or  guide  —  elastic 
fields  laterally.  The  basic  premise  in  designing  "surface  waveguides" 
is  that  a  strip  of  the  surface  on  which  the  characteristic  surface  wave 
velocity  is  lower  than  that  of  the  surroundings  will  tend  to  guide  sur¬ 
face  waves  along  its  length.  This  is  sometimes  implemented  by 
depositing  a  thin  layer  of  a  heavy,  soft  material  on  the  surface  in  the 
form  of  a  strip.  The  thin  layer  "loads"  the  substrate,  creating  a 
zone  in  which  the  surface  wave  velocity  is  lower  than  that  in  the 
unplated  region.  Conversely,  on  a  surface  fully  plated  with  a  light, 
stiff  material  except  for  a  strip  that  is  left  unplated,  the  strip  again 
acts  as  a  guide  since  the  surface  wave  velocity  is  higher  in  the  coated 
region  f^]1. 

Straight  surface  waveguides  of  the  type  described  above  have 
been  analyzed  by  Tiersten  [3],  who  represents  the  effect  of  the  sur¬ 
face  layer  by  a  1  mndary  condition  at  the  surface  of  an  isotropic 
substrate.  He  derives  from  his  model  an  equation  for  the  velocity  of 
the  "straight-crested"  (i.  e.,  plane  strain)  surface  waves  character¬ 
istic  of  the  plated  region.  He  then  appeals  to  a  concept  introduced  by 
Knowles  [4]  to  demonstrate  the  existence  —  in  the  same  plated 
region  —  of  a  more  general  class  of  surface  waves  referred  to  as 
"variable-crested"  surface  waves.  The  full  field  in  the  substrate  is 
synthesized  approximately  from  the  variable-crested  surface  waves 
characteristic  of  the  plated  and  unplated  regions,  by  applying  an 
approximate  matching  condition  at  the  boundary  between  the  regions, 

Reference  [2]  provides  a  fuller  discussion  of  these,  as  well  as  other 
techniques  for  guiding  surface  waves. 
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and  by  requiring  the  surface  waves  outside  the  strip  to  decay  exponen¬ 
tially  along  the  surface  away  from  the  strip.  This  leads  to  an  approxi¬ 
mate  dispersion  relation  for  the  guide,  i.  e.  ,  an  equation  relating  the 
frequency  and  velocity  of  disturbances  propagating  along  the  guide. 
Tiersten  and  Davis  [5]  have  extended  this  analysis  to  the  case  of  a 
curved  surface  waveguide. 

Freund  [6]  has  treated  a  structure  similar  to  that  treated  in 
[3]  using  a  more  exact  approach.  The  surface  of  the  guiding  strip 
considered  in  [6]  is  assumed  to  be  traction-free,  while  the  surface 
outside  the  guide  is  assumed  to  be  free  of  shear  stress  and  to  suffer 
no  vertical  displacement.  The  scattering  of  a  single  straight-crested 
surface  wave  incident  on  the  edge  of  a  semi-infinite  guide  at  arbi¬ 
trary  angle  is  solved  exactly.  An  approximation  to  the  dispersion 
relation  for  a  guide  of  finite  width  is  obtained  by  superposing  incident 
and  reflected  straight-crested  surface  waves  and  neglecting  at  one 
edge  of  the  guide  the  body  waves  generated  at  the  other  edge.  Fossum 
and  Freund  [7]  have  obtained  and  solved  the  dispersion  equation  for 
the  guide  by  an  independent  method  based  on  an  integral  equation,  and 
thereby  have  validated  Freund's  original  approach. 

Apart  from  the  problem  solved  by  Freund  in  [6]  as  an  inter¬ 
mediate  step  in  his  analysis,  there  are  no  known  exact  solutions  for- 
scattering  of  an  elastic  surface  wave  by  the  edge  of  a  surface  layer. 
Approximate  solutions  employing  juxtaposition  of  straight-crested 
surface  wave  modes  have  beer  obtained  by  Li  et  al  [8]  and  McGarr 


and  Alsop  [9]. 


In  contrast,  exact  solutions  do  exist  for  the  scattering  of  cer¬ 
tain  types  of  electromagnetic  surface  waves.  Of  particular  interest 
is  the  wot k  of  Kay  [10],  who  considers  the  scattering  of  an  electro¬ 
magnetic  surface  wave  by  a  discontinuity  in  a  property  of  the  surface 
known  as  normal  impedance.  An  exact  solution  is  obtained  by  using 
the  Wiener-Hopf  technique.  The  problem  is  closely  analogous  to  the 
present  one,  and  in  fact  has  served  as  a  "proving  ground"  for  some  of 
the  techniques  employed  here. 

To  complete  this  brief  summary  of  related  literature,  note 
should  be  taken  of  the  study  of  Koiter  [11]  of  the  static  problem  of 
load  transfer  from  a  stringer  into  a  plate.  The  problem  is  one  of 
generalized  plane  stress,  but  is  essentially  identical  to  that  of  a 
half- space  in  plane  strain,  half  of  whose  surface  is  plated  with  a  thin 
layer,  with  a  static  line  load  acting  on  the  edge  of  the  layer  and  paral¬ 
lel  to  the  surface.  The  effect  on  the  substrate  of  the  stiffness  of  the 
layer  is  modeled  in  precisely  the  same  manner  as  in  the  present 
analysis. 

Although  the  present  consideration  of  the  scattering  of  a  sur¬ 
face  wave  by  a  thin  layer  was  motivated  by  the  signal-processing 
applications  discussed  above,  the  results  of  the  analysis  may  also  be 
of  interest  in  connection  with  seismology. 

This  study  begins  in  Section  1  with  a  review  of  pure  surface 
waves  in  half-spaces  with  the  surface  either  fully  traction-free  or 
fully  covered  with  a  thin  layer.  In  Section  2  the  central  problem  is 
formulated  -  that  for  the  scattering  of  a  Rayleigh  wave  normally 
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incident  on  the  edge  of  a  semi-infinite  layer.  In  Section  3  the  problem 
is  solved  exactly  for  a  thin  massless  layer,  and  in  Section  4  it  is 
solved  by  a  perturbation  process  for  a  very  thin  layer  possessing  both 
mass  and  stiffness. 
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1  •  Pure  Surface  Waves 

For  the  purposes  of  this  study,  a  pure  surface  wave  is  a  solu¬ 
tion  of  the  elastodynamic  equations  of  motion  in  a  half- space  which 
propagates  in  a  direction,  parallel  to  the  surface  and  decays  exponen¬ 
tially  with  depth.  Such  waves  can  exist  in  a  half  space  with  the  sur¬ 
face  either  fully  traction-free  or  fully  covered  by  a  layer  of  a  differ¬ 
ent  elastic  material.  In  this  section  the  properties  of  pure  surface 
waves  are  reviewed,  with  emphasis  upon  a  method  of  approximating 
the  effect  of  a  surface  layer  in  the  case  where  the  layer  is  thin.  It  is 
necessary  to  understand  these  waves  before  the  problem  correspond¬ 
ing  to  the  scattering  of  a  surface  wave  can  be  formulated. 

Consider  a  half-space  occupied  by  an  elastic,  isotropic,  homo¬ 
geneous  solid  with  Lame  parameters  X  and  |jl  and  density  p,  covered 
by  a  layer  of  thickness  2h'  of  a  different  material  with  Lame  param¬ 
eters  X  and  p.  and  density  p'.  A  coordinate  system  is  established 
as  shown  in  Fig.  1,  with  a  y-axis  running  out  of  the  plane  of  the 
paper.  The  displacements  in  the  x-,  y-,  and  z-directions  in  the 
half- space  are  denoted  by  U,  V,  and  W  respectively,  and  in  their 
fullest  generality  are  functions  of  x,  y,  z,  and  time  t.  The  case  of 
plane  strain  is  assumed,  so  that  all  displacements  are  independent  of 
y,  and  VsO.  Further,  the  time  dependence  is  assumed  to  be  har¬ 
monic  with 

U(x,  z,  t)  =  u(x,  z)e-iu,t  , 

-id)  t 


W(x,  ■/.,  t)  =  w(x,  z)e 


,  oj  >  0  . 


(1.  1) 
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For  convenience,  the  amplitudes  u,  w  will  henceforth  be  referred  to 
as  displacements.  Since  actual,  physical  displacements  must  be 

t 

real,  those  corresponding  to  the  present  problem  are  either  the  real 
or  imaginary  parts  of  the  quantities  appearing  in  (1.  1). 

Under  the  foregoing  assumptions  the  elastodynamic  equations 
of  motion  will  be  satisfied  if  the  displacements  are  related  to  the  two- 
dimensional  Lame  potentials  $(x,  z),  Y(x,  z)  by^ 


u(x,  z) 
w(x,  z) 


_9£  91 

9x  9z  ’ 

-  9±  .  91 

9z  9x  ' 


(1.2) 


and  if  the  potentials  satisfy  the  reduced  wave  equations 

V2§+k^$  =  0  ,  ^ 


where 


V2Y+k2Y  =  0  ,  -oo<x<oo,  zsO, 
s 

2  2 

v2-0l+^i 

9x^  &z£ 


kd  =  Uj/cd'  ks  =  Uj/cs  ' 

c  ,  =  V(^+2[i)/p  ,  c  =  rJuTp  .  J 


(1.3) 


The  potentials  §  and  Y  are  referred  to  respectively  as  the  dilatation 
and  shear  potentials;  the  speeds  and  cg  are  respectively  the  dila¬ 
tation-  and  shear-wave  speeds. 

2 

The  stresses  are  related  to  the  potentials  by  the  formulas 

^ee  [12] 

2 

6For  convenience,  the  quantities  appearing  in  (1.  4)  are  referred  to  as 
stresses;  the  actual,  physical  stresses  are  either  the  real  or  imagi¬ 
nary  parts  of  these  quantities  multiplied  by  exp[-iU)t], 
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axx:';  ^[(2kd-ks  +  2^)$-2l^11 


'N 


a  =  u[2|^~ -(kZ+2-^);F], 

xz  L  oxoz  s  0z<2 


(1.4) 


a  =  uC(2k2-k2+2-^)H2-^:]  . 

zz  L  a  s  gz2  oxo  z 

The  pure  surface  waves  of  interest  here  are  solutions  of  (1.  3) 
which  decay  exponentially  with  z  and  satisfy  appropriate  conditions  at 
z  =  0.  Such  solutions  take  the  form 


whe  re 


i(x,2)  =  Ae-nd(k,z±ikx, 

T(x,Z)  =Ape-nS(k)z±ito!, 

n  (k)  =  yk2-kp,  k>ly  P  =  d,  s, 


(1.  5) 


and  where  A,p,  and  the  wave  number  k  are  constants  yet  to  be  deter¬ 
mined. 


Now  assume  that  the  layer  covering  t?:o  half- space  is  thin 
enough  that  the  displacements  within  the  layer  can  be  taken  approxi¬ 
mately  as  independent  of  z.  By  virtue  of  the  requirement  that  dis¬ 
placements  be  continous  across  the  interface  between  the  layer  and 
the  half- space,  the  displacements  in  the  layer  must  then  coincide  with 
those  in  the  ha’f-space  at  z  =  0. 

Consider  the  forces  acting  on  a  differential  element  of  the  layer 
of  unit  width  in  the  v-direction  and  of  length  dx,  as  shown  in  Fig.  2. 
Upon  neglecting  all  internal  stress  resultants  in  the  layer  except  those 
due  to  simple  stretching,  the  balance  of  momentum  in  the  x-  and  z- 
directions  requires 


-10- 
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Figure  2.  Differential  Element  of 
The  Surface  Layer  . 
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doi  2 

2h/-^-  +  0^m  (x,  0)  +  2h#p1l)  u(x,  0)=0, 


xz 


a  (x,  0)  +  2h'puAv(x,  0)  =  0  » 
zz 


(1.6) 


where  a'  is  the  axial  stress  in  the  layer.  The  stresses  ^xz(x,  0)  and 

a  (x  0)  are  those  exerted  by  the  substrate  on  the  layer.  The  axial 
zz' 

stress  a'  is  related  to  the  displacement  u(x,  0)  by 

,.  .  .  /.  X'+u' „  9u(x,  0) 

°<<x)=4^(TW)  3x  ' 

Substitution  of  (1.  7)  into  (1.6)  yields 

s2 


(1.7) 


o  (x,  0)  +pN  9  +  pMu|x,  0)  =  0 , 1 

XZ 


9x 

0  (x,  0)  +  jjMw(x,  0)  =  0  , 

z  z 


(1.8) 


..  8h'u'.  X'+u\  w_2h^L 


whei*e 

N  -- 

U  'A'+Zu'"  —  M  J 

Note  that  N  and  M  are  measures  of  the  layer's  stiffness  and  mass. 


respectively. 

Equations  (1.8)  may  be  thought  of  as  boundary  conditions  to  be 
satisfied  by  the  fields  in  tie  substrate.  Note  that  in  the  absence  of  a 
covering  layer,  M=N  =  0  and  (1.8)  reduce  to  the  boundary  conditions 
appropriate  to  a  traction-free  surface.  Tiersten  [3]  deduced  (1.8) 
from  a  more  general  development  in  which  Mindlin's  plate  equations 
[13]  were  taken  as  a  starting  point^.  He  argued  that  for  thin  plates 


iThis  relationship  is  derived  from  the  three-dimensional  stress- strain 
law  under  the  assumption  of  plane  strain  in  the  y-direction  and  vanish¬ 
ing  of  normal  stress  in  the  z-direction. 

2Achenbach  and  Keshava  [14]  used  a  similar  approach  which  included 
fhe  effects  of  axial  stretching,  bending,  shear,  and  rotatory  inertia  in 

the  layer. 
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all  terms  except  those  involving  axial  stretching  and  inertial  reaction 
are  negligible. 

Substitution  of  the  potentials  (1.  5)  into  (1.  2)  and  (1.  4),  and 
thence  into  the  boundary  conditions  (1.  8),  yields  a  pair  of  homoge¬ 
neous.-  linear,  algebraic  equations  in  A  and  A|3.  The  condition  for 
the  existence  of  a  non-trivial  solution  is  that  the  determinant  of  the 
coefficient  matrix  vanish.  This  is 

R(k)  -  Nk2k^ns(k)  +  M{k2[ng(k)  +  nd(k)] 

+  Nk2[k2-  ns(k)nd(k)]HM2{ns(k)nd(k)  -  k2]  =  0  (1.  9} 

where 

R(k)  =  (2k2-  k2)2-  4k2ng(k)nd(k)  .  (1.10) 

Each  value  of  k  satisfying  (1.  9)  determines  an  elastic  surface 
wave  through  (1.5),  provided  k>kg  (note  that  kg>kd,  provided  X>0 
and  u>0  —  see  (1.  3))  and  provided  (3  has  the  now  common  value  re¬ 
quired  by  either  of  the  boundary  conditions  (1.8).  If  the  time  depen¬ 
dence  is  restored  in  the  potentials  (1.  5),  their  x-  and  t-dependence, 
and  hence  that  of  all  field  quantities,  is  of  the  form 

e±ikx-iujt=  e±ik(xTct)  {1 

where  c-tu/k.  Thus  the  surface  wave  (1.5),  while  decaying  exponen¬ 
tially  with  z,  propagates  unchanged  in  form  with  speed  c  in  the 
positive  or  negative  x-direction  according  to  whether  the  positive  or 
negative  sign  is  used  in  (1.  11). 

If  the  surface  layer  thickness  2h'  vanishes,  M=N=0,  the  sur¬ 
face  becomes  traction-free,  and  (1.9)  reduces  to 
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R(k)  =  0  .  (1.  12) 

Equation  (1.  12)  determines  the  propagation  speed  of  Rayleigh  waves, 
an  extensive  discussion  of  which  may  be  found  on  pp.  307  -  309  of 
Love  s  text  [15],  Achenbach,  on  p.  190  of  [16],  proves  that  there  is 
only  one  root  of  (1.  12)  greater  than  kg.  This  root  is  denoted  by  kR 
and  referred  to  as  the  Rayleigh  root  or  the  Rayleigh  wavenumber.  If 
(1.  12)  is  multiplied  by  (U  the  resulting  equation  for  cD  =  ou/k-,  and 
hence  cR  itself,  is  independent  of  U).  The  Rayleigh  wave  is  thus  said 
to  be  non-dispersive,  i.  e.  ,  its  propagation  speed  cR  is  independent  of 
frequency  CU. 

In  the  more  general  case  Eq,  (1,9)  rnay  have  one  or  two 
roots  [3],  The  smaller  one,  which  will  be  denoted  kT  and  referred 
to  as  the  extended  Rayleigh  root,  reduces  continuously  to  kD  as  M  and 
N  go  to  zero,  while  the  other  root  goes  to  infinity  in  this  limit.  If 
N  /  0  or  M  ^  0,  the  speed  corresponding  to  the  extended  Rayleigh  root 
is  frequency-dependent,  so  the  presence  of  the  layer  results  in  sur¬ 
face  waves  which  are  dispersive. 

A  special  case  of  interest  is  that  of  a  massless  layer,  for  which, 
bV  (1-8),  M  =  0,  and  (1.9)  reduces  to 

Rn'(k)  =R(k)-Nk2k2n  (k)  =  0  .  (1.13) 

s  s 

It  is  shown  in  Appendix  A  that  this  equation  also  has  but  one  admissi¬ 
ble  root.  This  root  is  less  than  k„ ,  reduces  continuously  to  k_,  as 
N-»0r  and  corresponds  to  a  dispersive  surface  wave  which  propagates 
faster  than  a  Rayleigh  wave  but  slower  than  a  shear  wave.  This  is 
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in  keeping  with  the  intuitive  expectation  that  a  stiffening  layer  would 
increase  the  characteristic  surface  wave  speed. 

It  must  be  noted  that  if  the  sole  objective  of  the  present  paper 
were  a  study  of  pure  surface  waves,  it  would  not  be  necessary  to  make 
any  simplifying,  approximative  assumptions  regarding  the  displace¬ 
ments  in  the  layer.  Rather,  the  layer  could  be  treated  as  a  second 
elastic  continuum,  and  exact  solutions  to  the  field  equations  in  both 
the  layer  and  the  substrate  could  be  obtained,  subject  to  the  appropri¬ 
ate  continuity  conditions  at  the  interface,  boundary  conditions  at  the 
traction-free  surface,  and  decay  conditions  in  the  substrate.  Auld 
summarizes  this  topic  on  pp.  97  -  103  of  [17];  Ewing,  Jardetzky,  and 
Press  discuss  it  on  pp.  189  -  223  of  [  18].  Tier sten  [3]  compares  the 
solutions  of  the  exact  dispersion  equation  with  those  of  (1.  9),  conclud¬ 
ing  that  the  results  agree  extremely  well  for  layers  less  than  about 
one  sixth  of  a  wavelength  thick,  and  reasonably  well  up  to  about  half  a 
wavelength.  The  reason  the  exact  analysis  is  not  used  here  is  that  it 
does  not  allow  the  effect  of  the  layer  to  be  fully  specified  by  a  simple 
pair  of  boundary  conditions  such  as  (1.  8),  The  treatment  of  the  effects 
of  the  layer  through  simple  boundary  conditions  is  essential  to  the 
tractability  of  the  scattering  problem  which  is  the  main  objective  of 


this  work. 


2.  Formulation  of  the  Problem 


An  elastic,  isotropic,  homogeneous  solid  occupies  the  half¬ 
space  -oo<x,y<oo,  z&O  .  A  thin  layer  covers  and  is  bonded  to  the 
surface  z  =  0,  -oo<y<oo,  x>0,  as  shown  in  Fig.  3.  Notations  for 
coordinates,  material  properties,  and  field  quantities  are  exactly  as 
in  the  previous  section.  The  assumptions  of  plane  strain  and  harmonic 
time  dependence  of  the  form  e  are  maintained,  so  that  the  problem 
may  be  formulated  in  terms  of  the  time-independent  Lame  potentials 
§  (x,  z),  Y(x,  z).  The  reduced  wave  equations  governing  these  poten¬ 
tials  are  (1.3),  the  displacement-potential  relations  are  (1.2),  and 
the  stress-potential  relations  are  (1.4).  The  boundary  conditions  for 
the  unplated  surface  are 


(x,U)  =  0  , 


(2.  1) 


a  <x,  0)  =  0  ,  -co<x<0  , 

XZ  ] 

and  for  the  plated  surface,  in  accordance  with  the  assumptions  and 
results  of  the  previous  section. 


a  (x,  0)  +  |iMw(x,  0)  =  0  , 

zz 


a  (x,  0)  +  pN 


xz 


8  u(x,  0) 
8x2 


+  pMu(x,  0)  =  0  ,  0<x<oo  . 


(2.  2) 


The  potentials  §,¥  are  required  to  be  twice -continuously 
differentiable  in  the  interior  of  the  half- space.  The  quantities  appear¬ 
ing  in  the  boundary  conditions  (2.  1),  (2.  2)  are  required  to  be  continu¬ 
ous  up  to  the  respectively  appropriate  segment  of  the  boundary,  ex¬ 
cepting  the  origin  x  =  z  =  0  .  Certain  further  conditions,  which  will  be 
explicitly  specified  in  the  next  section,  are  required  to  be  satisfied  by 
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Figure  3.  Coordinates  for  a  Half-Space  Covered 
by  a  Semi-Infinite  Thin  Layer. 
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various  field  quantities  as  the  origin  is  approached,  these  conditions 
having  the  effect  of  ruling  out  point  loads  and  higher-order  singulari¬ 
ties  at  the  origin. 

The  potentials  $,  ¥  are  required  to  be  of  the  form 

$  (x,  z)  =  A0cpR+(x,  z)H(-x)  +  ArcpR_(x,  z)H(-x) 

+  A  cp™  ,  (x,  z)H(x) +cp(x,  z)  , 

t  T+  \  (2.  3) 

■  tx,  z)  =  AQt|/R+(x,  z)H(-x)  +  Ar1'R_(x,  z)H(-x) 

+  At^T+^X’  +  ^x»  z)  •  > 


In  Eq.  (2.3),  cpR± ,  fR± ,  cpT+ ,  and  i|>T+  are  normalized  potentials 
representing  positive-  and  negative -traveling  Rayleigh  waves  and  the 
positive-traveling  extended  Rayleigh  wave  corresponding  to  the  layer, 


and  are  given  by 

°PR± 


where 


=  e-nd<kR)z±2kRx,  *OJ.=  #De-n3(kR)a:tikRX, 


A 


cpT+ 


R±  R 

=  e-nd(kT*a  +  iklx,  *T+=  pTa-”S(kT)z+ikTx, 


,.2  ,2 

2kR"ks 

iR  2ikRns(kR)  ’ 


(2.4) 


2kT'k3*|iMnd^kT^ 
”  2ikTng(kT)-i|iMkT 


J 


The  constant  Aq  is  given,  the  constants  Af  and  Afc  are  unknown, 

H(x)  is  the  unit  step  function,  and  cp  and  f  are  required  to  satisfy  the 


radiation  conditions 


ik  ,r 

cp(x,z)~F(e)^r-  , 

ik  r 

iji(x,  z)~G(e)e  n  , 


(2.  5) 
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/  2  2  -1 

as  r-*oo,  uniformly  in  0  for  Os0£ir  .  Here  r  =  Vx  ,  6  =  Tan  (z/x), 
F(0)  and  G(0)  are  unknown  functions  of  0,  and  the  symbol  n~" 
denotes  asymptotic  equality  in  the  usual  sense. 

The  motivation  behind  (2.  4)  will  be  discussed  presently,  but 
first  a  matrix  notation  will  be  introduced  which  not  only  simplifies 
some  of  the  subsequent  manipulations,  but  consolidates  the  fundamen¬ 
tal  concepts  of  the  analysis.  Twice -underlined  letters  will  denote 
2x2  matrices  of  partial  differential  operators,  and  underlined  letters 
will  denote  column  vectors  with  two  components,  the  first  of  which  is 
a  dilatational  potential,  and  the  second,  the  corresponding  shear 
potential.  Thus  the  reduced  wave  equations  (1  3)  become 


where 


L  = 


L  £  =  0 

^2+kd  0  ^ 

*(x,z) 

\ 

o  v*+k j 

T(x, z) 

J 


The  boundary  conditions  (2.  1)  and  (2.  ?.)  become 

II  =  _0  ,  z  =  0  ,  -ao  <x<0  ,  ^ 

(R  +  NS  +  MT)£  =  0  ,  z  =  0  ,  0<x<oo  , 


where 


R  = 


4k2-k2+2-^ 

d  s  9z2  oxoz 
a2 

9x3z 


f  0 


A 


»■ 


dJ _ 

\^3x^  Ox^dz J 


J 


(2.6) 


(2.7) 
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The  required  form  of  the  solution  (2.  3)  becomes 


where 


1  =  A0C£R+H{-x)+Ar^R^H(-x)  +  AtceT+H(x)+^ 


£ 


R± 


e~nd(kR)z 

±pRe'ns(kR)z 


e±lkRx  , 


)  (2.8) 


£t+  = 


fe-nd(kT,z 

PTe'ns(kT)z 


:ikTX, 


cp  = 


Cp(x,  z) 

z) 


and  where  tp(x,  z)  and  i|r(x,  z)  must  satisfy  the  radiation  conditions 
(2.5). 


Now  it  can  be  seen  that  in  the  assumed  form  (2.  8)  each  of  the 
first  three  terms  represents  a  surface  wave  propagating  in  either  the 
positive  or  negative  x-direction  and  existing  only  for  either  positive  or 
negative  x.  The  first  term  is  the  incident  wave  and  has  amplitude 
Aq  .  The  remaining  terms  taken  together  constitute  the  scattered 
wave  and  propagate  outward  from  the  origin.  These  terms  are 
referred  to  respectively  as  the  reflected  wave,  the  transmitted  wave, 
and  the  radiated  wave;  the  ratios  A^Aq  and  A^/Aq  are  referred  to 
respectively  as  the  reflection  and  transmission  coefficients. 

Note  that  the  radiated  wave  must  contain  discontinuities  across 
x  =  0  which  cancel  those  due  to  the  surface  wave  terms  of  (2.  7),  so  as 
to  render  the  full  field  twice-continuously  differentiable.  Morgan 
and  Karp  [19]  have  solved  an  electromagnetic  problem  involving  sur¬ 
face  waves  by  using  a  similar  representation,  in  which  the  full  field 
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is  separated  into  various  parts,  some  representing  surface  waves, 
one  representing  a  radiated  wave,  and  all  possessing  discontinuities. 
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3 .  Solution  for  a  Massless  Layer 

It  is  now  assumed  that  the  layer's  density  p'  vanishes. 
According  to  (1.  8),  this  causes  M  to  vanish.  Under  this  condition  the 
governing  differential  equations  (2.  6)  remain  the  same,  but  the  bound¬ 
ary  conditions  (2.  7)  reduce  to 


R  <1  =  0  ,  z  =  0  ,  -oo<x<0  ,  > 

=  ~  ”  l  (3.1) 

(R+NS)i=0,  z  =  0  ,  0<x<co  J 

The  assumed  form  (2.  8)  and  the  radiation  conditions  (2.  5)  remain,  but 
the  wave  number  krj,  appearing  in  (2.  8)  is  now  a  root  of  the  simplified 
form  (1.  13)  of  the  more  general  dispersion  relation  (1.  9). 

The  order  conditions  at  the  origin  are  taken  as 

a  (x,0)=O(l),  > 

X  A 

3a  (x, 0)  , 

— Xgx -  =  0(x  ^  )  as  x-*0  ,  pi  <1  ,  )  (3.2) 

2 

a  (x,  0)  +uN  -  ^  =  0(  jx|"^s)  as  x-*0  ,  pa<l  . 

XZ  3x^  J 

The  solution  to  the  problem  specified  by  (2.  6),  (3.  1),  (2.  8), 

(2.  5)  and  (3.  2)  will  be  obtained  in  closed  form  in  this  section. 

A 

It  is  convenient  to  define  a  new  function  £(x,  z)  by 

£(x,  z)  =  A0C£R+(x,  z)  +£(x,  z)  .  (3.3) 


Recall  that  is  the  amplitude  of  the  incident  Rayleigh  wave.  Thus 

A 

for  x<0,  £  comprises  the  scattered  wave,  while  for  x>  0  it  consists 
of  the  scattered  wave  minus  a  positive -traveling  Rayleigh  wave  of 

A 

amplitude  AQ.  In  either  case,  JP  consists  only  of  outwardly- propagat¬ 


ing  waves. 
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The  differential  equations  governing  $  are  obtained  by  substi¬ 
tuting  (3.  3)  into  (2.  6),  and  by  noting  that,  according  to  (1.  3}  and  (2.  6), 
jL<£R+-0.  This  yields 

L  cp  =  0  ,  - oo <x <oo  ,  zsO  .  4) 


The  boundary  conditions  for  £  are  obtained  by  substituting  (3.3)  into 
(3.  1),  and  by  noting  that,  according  to  (1.  8),  (1.  2),  (1.  4)  and  (2.  7), 


R  cpR+  =  0.  This  yields 


R  £  =  0  ,  z  =  0  ,  -oo<x<oo 


ArtNik  k 
0  s 


£  j°|elkRX, 


z  =  0  ,  0<x<oo  , 


(3.  5) 


(R  +  N3)£=  2 

where  use  has  been  made  of  the  definition  of  S  in  (2.  7)  and  that  of 


PR  in  (2.  4). 

This  problem  will  be  solved  with  the  aid  of  a  Fourier  trans¬ 
form  on  x.  Let  f(x)  be  a  complex-valued,  absolutely  integrable  func¬ 
tion  of  the  real  variable  x  which  satisfies 

|f(x)|<Ke'5K  |x|>xo>  (3'6’ 

for  some  positive  constants  K,  6  and  Xfl.  The  exponential  Fourier 
transform  of  f(x)  is  defined  by 

CD 

f'(H)  =  Jf(x)e"WXdx  ,  -6<Im{K}<6,  (3-7) 

-CO 


and  the  inversion  formula  is 

oo 

f(x)  =  Jf  (K)c1KXdK  ,  -00<X<CD  . 
-C X> 


(3.8) 


When  considered  as  a  function  of  the  complex  variable  x,  the 
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transform  f  (h)  is  analytic  in  the  domain  -6<Im[n}<6  (see  [20], 
p.  338). 

It  is  clear  from  (2.  8)  that  for  fixed  z  the  potentials  §,  ¥  will 
not  satisfy  the  decay  condition  (3.  6)  for  any  as  long  as  the  wave 

numbers  k^,  k^,,  kg  and  k^  are  real.  However,  application  of  the 
Wiener-Hopf  technique  (which  will  be  explained  presently)  requires 
that  the  transforms  be  analytic  in  strips  of  the  complex  K-plane  such 
as  that  defined  in  (3.  7).  The  difficulty  is  resolved  by  temporarily 
assuming  that  the  frequency  uu  has  a  small,  positive,  imaginary  part, 
while  all  propagation  speeds  remain  real.  This  is  a  standard  artifice; 
it  is  discussed,  for  example,  on  p.  28  of  [21].  After  the  solution  is 
obtained,  the  imaginary  pare  of  uu  will  be  taken  to  vanish. 

Now  since  (£  13  required  to  contain  only  outwardly  propagating 
waves,  the  assumption  that  lm{(ju}>0  (and  consequently  that 
Irn  (k^l  >  0  ,  tv  =  d,  s,  T,  R)  renders  $  exponentially  decaying  in  |x|  as 
required  by  (3.6),  with  decay  constant  6  satisfying  0<6 <Im ( k^} 

(cf.  (3.3),  (2.8),  (2.5)). 


Formal  application  of  the  transform  (3.  7)  to  the  differential 
equations  (2.  6)  yields 


■\ 


.2* 


-  „*<»)$<>..  z)  =  0  , 
dz*  d 


a2 T(h  Z)  2  *  >  (3- 10> 

-  7^  -  n  (h)$(h,  z)  =  0  ,  -6<Im{H}<6, 

dt  s 


where 


J  2  .  2'  J  2  .2' 

n,  =  VK  -  k  ,  ,  p  =v*-k 
a  a  s  s 


J 


) 
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The  functions  ng(x)  and  n^fx)  are  rendered  single- valued  in  the  entire 
H-plane  and  analytic  in  the  strip  -6<Im{x}<6  by  taking  branch  cuts 
as  shown  in  Fig.  4,  and  by  requiring  n  (x)  and  n  ,(x)  to  approach 
+  oo  as  x-oo  along  the  positive  real  x-axis.  The  solutions  of  (3.  10) 
are  then 

$(x,z)  =  A(x)e"nd(*,z+  A'(x)end(x)z , 


?(x,z)  =  B(x)e‘ns(H)z+  B'(x)eVK,*f  -6<Im{x]<6,  j 


(3.  11) 


where  A(x),  A'(x),  B(x),  and  B'(x)  are  unknown  functions  of  x.  How¬ 
ever,  the  coefficients  A'(x)  and  B'(x)  must  vanish;  otherwise  the 
transforms  would  grow  exponentially  with  |x|  for  z>0,  and  the  inver¬ 
sion  integrals  (3.8)  would  not  converge.  This  leaves 


&X.Z)  =A(x)e-nd(H)2f 
£{x,z)  s  B(x)e‘ns(H)z,  -6<Im{x)<ft 


(3.  12) 


The  boundary  condition  (3.  5),  when  written  in  component  form, 
becomes 

2  2* 

(2kd‘kf+2Tr)^+  ZlkJz  =  °*  z=0»  (3.13a) 


2*  2 

2  9xDz  "  \  s  +  2  2  =  z  =  °»  -oo<x<0  ,  (3.13b) 

8z 


(2  K5?  +  •(k!+27-2tN7X-)* 


>2 

dx 


dz  dx  dz‘ 


-  i[^0^'kskRe^^X  '  *  =  ®  *  0<x<oo 


(3.  13c) 


Note  that  (3.  13a),  which  refers  to  the  normal  stress  at  the  surface, 
holdsforall  values  of  x  at  z  =  0.  This  circumstance  arises  from  the 


Figure  4..  Complex  h- Plane. 
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fact  that  the  mass  of  the  layer  is  neglected  at  present,  so  that  M  =  0  in 
(2.  2).  In  contrast,  the  conditions  (3.  13b,  c),  which  relate  to  the 
normal  tractions  at  z=0,  take  different  forms  according  as  x<0  or 
x>0.  The  fact  that  only  one  of  the  boundary  conditions  changes  form 
at  x  =  0  makes  possible  the  application  of  the  Wiener-Hopf  procedure 
in  its  simplest  form. 

It  is  now  convenient  to  define  functions  f(x)  and  g(x)  on  (0,oo) 
and  (-oo,0),  respectively,  by  setting 


2* 

f(x)  =  z=0>  O<x<00’ 

az 


.2 


g(x)  =  (; 


A2  A3 


k2+2-iU:+N— %- 

S  Q  “  r\  "  r\ 

o z  dx  dz 


')$ , z=0  ,  - 


oo<x<0 . 


(3.  14) 


Apart  from  a  factor  fi,  f(x)  and  g(x)  are  the  unknown  values  of  a  and 

X  z 

2  2 

a  +  jjN8  u/3x  under  the  layer  z  =  0,  0<x<oo  ,  and  on  the  free  sur- 
xz 

face  -oo<x<0  ,  z  =  0  ,  respectively. 

Formal  application  of  the  transform  (3.  7)  to  (3.  13)  and  (3.  14), 
followed  by  use  of  (3.  12),  yields 


(2k2-  k2)A(n)  -  2iH.ng(H.)B(n)  =  0  , 

-2i?ui  (k)A(h)  -  (2k2-  kf)B(K)  =F‘(h)  , 
s  s 

[-2iKn  (k)  -  iNH'']A(H)  -  L  2h2- k2+NK2n  (k)]B(k) 
s  s  s 


> 


AnNk2kD 
0  a  R 

2(x-kR) 


+  G+(h),  -6<Im{K}<6  , 


(3.  15) 


J 
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where 


oo 

F~(k)  =  Jf(x)e  iHXdx 
0 

0 

G+(h)  =  Jg(x)e‘1Hxdx 

-oo 


■N 


J 


(3.  16) 


The  functions  f(x)  and  g(x)  as  defined  by  {3.  14)  are  exponen¬ 
tially  decaying  for  large  positive  and  negative  x,  repsectively,  with 
decay  constant  6.  Moreover,  since  F  (h)  and  G+(h)  are  half- range 
transforms  of  f(x)  and  g(x),  F"(h.)  and  G+(h)  are  analytic  functions 
of  h  on  the  half- planes  Imfx}<6  and  Im{x]>-6,  respectively  (see 
p.  12  of  [21]).  The  superscripts  on  F  (x)  and  G+(x)  are  thus  seen 
to  be  suggestive  of  the  half-planes  on  which  the  functions  are  respec¬ 
tively  analytic. 

After  A(x)  and  B(h)  are  eliminated  from  (3.  15)  there  remains 
the  single  equation 


„ ,  .  _  ,  A„Nk  k„ 

F‘(H)  =  fw[G  M+ 

where 

R(x)  =  (2x2- k2)2- 4x2n  (x)n  ,(x)  , 

3  S  Q 

R'(x)  =  R(x)  -  Nk2x2n  (x)  . 

s  s 


(3.  17a) 


(3.  17b) 


It  is  important  to  note  that  the  domain  of  validity  of  (3.  17a)  is  the 
strip  -6<Lm{x}<<5. 

Now  define  D(h)  by 

2(kZ-  k2)R* (x) 

D(h)  =  - I - - - ,  -6<Im{x}<6  . 

Nk  n  (x)R(x) 
s  s 


(3.  18) 
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It  is  shown  in  Appendix  B  that  D(h)  admits  the  representation 


where 


d(k)  = 

D+(h) 

)  (3.19) 

CDTifi 

D^h)  =  exp{- ~rj  L9-gI:^L3da}  ,  -6<Im{Hl<6 
-crisis  y 


and  that  D^(x)  and  D  (x),  by  trivial  analytic  continuations,  have  the 
same  domains  of  analyticity  as  G+(x)  and  F  (x),  respectively.  In 
view  of  (3.  18)  and  (3.  19),  Eq.  (3.  17a)  may  be  rewritten  as 


NkV(H)D‘(*)F'(>0  AQN)ARDt(kR) 


2<lV  kd» 


2(x-kR)ng+(kR) 


1-.+  /  i  A-Nk^k^  %  D+(kr,)-1 

°  p  WO.  <al  t  0  »  f.rsjo  - -A],  -6<lmtk]<6, 

n+<*>  2(*-kR)  LVW  ns+(kR>- 

s 


(3.  20) 


where  n^x)  =  Jn±k  '  .  The  functions  n^x)  are  rendered  single- 
s  K.  s 

valued  in  the  entire  complex  x-plane  and  analytic  in  the  strip 
-*<Im(x}<6  by  taking  branch  cuts  as  shown  in  Fig.  4,  and  by  re- 

X  _  . 

quiring  both  n  (k)  and  n  (h)  to  approach  +co  as  h^oo  along  the  posi- 
s  s 

tive  real  axis.  Note  that  the  superscripts  have  the  same  meanings  as 
above. 


Equation  (3.  20)  is  said  to  be  of  Wiener-Hopf  type.  The  method 
used  to  solve  it  is  referred  to  as  the  Wiener  -Hopf  technique,  and  is 
outlined  on  p.  37  of  [21].  Define  functions  H±(h)  by 
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NkV(x)D"{H}F“(H)  AnNk\  D+(kR) 

H‘(h)  =  5 . 5- . -,--j - 0  8  R_,  R-  ,Im{xj<6, 


2"V  $ 


2 

~+,  ^+/  v  AnNk  kn  rn4,  .  D  (k 

n+M  =  D  -M  +  jknrfPV^ ~ 6 


2(K-kR)ns(kR) 


D+{k>. 


>(3. 


21: 


V*1 


The  function  H+(x),  by  a  trivial  analytic  continuation,  is  analytic  for 
Im{x}>-fi,  while  H~(x)  is  similarly  analytic  for  Jxn{xl<6.  Accord¬ 
ing  to  (3.  20},  there  is  a  common  strip  of  analyticity  in  which 


H+(x)  =  H”(x) ,  -6<im{x}<fi 


(3.  22) 


It  follows  that  there  is  an  entire  function  H(x)  such  that 


f  H" (h)  ,  Im{x}<fi  ,  1 

H(k)  =  ,  (3-23) 

Ih  (x) ,  Im{x}>  -6  .  J 

Now  it  will  be  proved  that  both  H  (x)  and  H+(x)  tend  to  zero  as 
|x|  -♦oo  in  the  appropriate  half-plane.  The  order  conditions  (3.  2)  imply 
the  existence  of  positive  constants  L  and  e  such  that 


oxb(x,°)|<l 


,  0<xse  , 


*\ 


do  (x,  0) 

-  |  <Lx"Pl  f  o<x££  , 


[axs(x,  0)  +  pN^%^  |  CL  \x\-te  ,  -e^x<0  . 


>  (3.24) 


3x 


J 


The  exponential  decay  resulting  from  the  complex  frequency,  together 
with  the  assumed  continuity  of  the  solution,  imply  the  existence  of  a 
constant  K>0  such  that 
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|axz(x,0)|<Ke  6x ,  e^x<oo, 

3ct  (x,  0) 

|  ^ - |<Ke  x,  esxCoo  ,  ) 

kx  (X»Q)  +  |-tN—  |<Ke^x  ,  -oo<x£-e. 

XZ  dx  ^ 


(3.  25) 


Substitution  of  (1,  2)  and  (1.  4)  into  (3.  14),  thence  the  result  into  (3.  16), 


yields 


oo 


|iF  (x)  =  fa 


*j  xz 
0 


(x,0)e"1HXdx, 


uG+{k)  =  J  [(J  (x,  0)  +^iN— -  ^  Q)le'1HXdx  . 

-ooL  XZ  3xZ  J 


(3.  26a) 

(3.  26b) 


One  integration  of  (3.  26a)  by  parts,  followed  by  use  of  (3.  24)  and 
(3.  25),  leads  to  the  following  asymptotic  estimates  for  F(x)  and 
G+(x): 


IhF“(k)|  = 


xz 


(x,  0) 


lXXl 


-IX 


GO 

0+ 


1  f  daxz(x’  0)  -i XX 

+  77  — ^ - e  1Kxdx 

xx  J  ox 

0 


5  Fl +  e[JLx~Pldx  +  JKe"6Xdx] 

1  1  I  I  0  e 

M3.  27) 

=  0(  |wj  )  as  |x|-*ao,  Im{x}<6, 

0  -e 

j|iG‘f(x)|  £  J  L|xj~^2dx  +  J  Ke^Xdx 
-  e  -oo 

=  0(1)  as  |x|-oo,  Im{x}>  -6  . 
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It  is  proved  in  Appendix  B  that 


|d+(k)|-i  as  |x,|~oo,  Im(K}>-6, 
|D~ (k)|  —  1  as  |h|  -*oo,  Im f k] <  5  . 


(3.  28) 


Use  of  (3.  27)  and  {3.  28)  in  (3.  21)  yields 

H  (k)-*0  as  |?i|-‘oo,  Im{K}<6  , 

4* 

H  (h)"*0  as  |k|  oo  P  Im  {  r  }  >  -  5  . 

Therefore  (3.  23)  and  a  trivial  extension  of  Liouville's  theorem  (p.  125 
of  [22]  )  imply  that 

H(k)  s  0  .  (3.29) 

In  view  of  (3.  29),  (3.23),  and  (3.  21),  F  (x.)  and  G+(x.)  are 


given  by 


F‘(h)  = 


AOkR^ks"kd^D  *kR^ 
n”  (h)  D~  (?0(H-kR)n^(kR) 


,  Im[ k ] <6  , 


AoNkfkR  r.  "sWD'(kR>i . I 


Gt(k)  =  - 


D  (H)ns(kR)- 


(3.  30) 


Substitution  from  (3.30)  into  (3.  15)  yields  three  equations  in  A(x)  and 
B(h).  Any  one  of  these  may  be  considered  as  redundant  and  the  re¬ 
maining  two  solved  for  A(h)  and  B(h)  to  give 


A(k)  - 


A0NkRksD+(kR)Kns(H)ns(K) 

in+(kR)(H-kR)R*(K)D+(x) 

2A0kRD+<kR)(ks-kd>Kns(H) 
in^(kR)n^(K)(H-kR)R(x)D  (h) 


t 


(3. 31a) 


(3.  31b) 


-32- 


B(k)  = 


A0NkRkV(kR)(2K2-k^) 

2n+(kR)(*-kR)R*(K)D+(x) 


(3.  31c) 


AOkRD  (kR^ks-kd,(2K2'k^ 

— -  ■, -6<Im{x}<6.  {3.  31d) 


ng(kR)n^{H){K-kR)R(H)D"(;t) 

The  alternate  fcurms  are  obtained  by  using  (3.  18)  and  (3.  19). 

Application  of  the  Fourier  inversion  formula  (3.  8)  to  (3.  12) 

yields 

CD 

cp(x,  z)  =  A(x)e"nd^z  +  1Kxdx  ,] 


-oo 


CO 


hx.z)  =  jB(*)e‘na(*,z  +  i*Itd* 


(3.  32) 


-co 


as  formal  candidates  for  a  solution  to  the  original  problem  for  cp, 
Several  steps  are  now  necessary  to  prove  that  the  functions 
cp,  $  appearing  in  (3.  32)  are  properly  defined  and  actually  satisfy 
(3.  4)  and  (3.  5).  These  steps  are  left  to  Appendix  C.  It  is  also  proved 
there  that  the  asymptotic  forms  of  F'(h)  and  G+(x)  as  |h|-cx>  are  in 
fact  as  specified  in  (3.  27). 

In  older  to  demonstrate  that  (3.  32)  contain  the  surface  waves 

lequired  by  (2.  8),  the  integrals  in  (3.  32)  must  be  regarded  as  complex 

contour  integrals  along  the  real  H-axis  in  the  complex  H-plane  (see 

Fig.  4).  The  singularities  in  the  integrands  in  (3.  32)  are  due  to  the 

various  singularities  of  D  (h)  and  D  (x),  the  branch  points  of  n  (x) 

s 

and  n^n),  and  the  simple  zeros  of  the  denominators  associated  with 
A M  and  B(x).  It  is  shown  in  Appendix  B  that  D+(x)  has  branch 
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points  at  -k  and  -kd,  no  zeros,  and  a  simple  pole  at  -kT>  while 

D"(k)  has  branch  points  at  kg  and  k^,  no  zeros,  and  a  simple  pole  at 

k  .  By  inspection  of  (3.  17b)  it  is  seen  that  R(h)  and  R  (x)  have 
R 

branch  points  at  ±kg  and  ±kd,  and  it  is  proved  in  Appendix  A  that 
their  only  zeros  are  located  at  =tkR  and  ±k^-,,  respectively.  Hence 
the  singularities  of  the  integrands  in  (3.3  2)  are  branch  points  at  ±kg 
and  ±kd>  and  simple  poles  at  ±kR  and  ±k,j,.  Ihese  are  shown  in 
Fig.  4.  That  they  are  not  on  the  real  axis  is  due  to  the  assumed  small, 
positive,  imaginary  part  of  uj.  The  branch  cuts  are  taken  as  shown 
in  Fig.  4.  The  segments  of  these  cuts  between  the  branch  points  and 
the  imaginary  axis  will  lie  along  the  real  axis  in  the  limit  as  Im{w]  -*0. 

Now  consider  an  integration  around  the  closed  contour  formed  by 
the  quarter-circular  arcs  of  radius  R  and  the  appropriate  segments 
of  the  branch-line  contour  and  the  real  axis,  as  shown  in  Fig.  4.  If 
the  integrals  over  the  quarter-circular  arcs  vanish  in  the  limit  as 
r-.qo,  then  by  Cauchy's  theorem  the  integral  over  the  real  axis  may 
be  replaced  by  that  over  the  branch-line  contour  plus  the  contribution 
from  the  residues  at  the  poles  k,j,  and  kR.  It  can  easily  be  shown 
that  this  is  in  fact  the  case  when  x>0.  The  residues  are  evaluated 
using  (3.  3  la,  c)  and  (3.  32).  With  the  aid  of  (3.  17b)  and  (3.  3)  it  is 
seen  that  the  contribution  due  to  the  pole  at  kR  is  equal  to  the  inci¬ 
dent  wave  A(^PR+(x,  z)  but  for  its  sign.  Thus  substitution  of  the 
result  of  the  contour  shift  into  (3.  3)  yields 


-34- 


where 


§  (x,  z)  =  AtcpT  +  {x,  z)  +  cp(x,  z)  , 

z)  =  At^T+(x,  z)  +  i|f(x,  z)  ,  x>0  , 


A 


t 


A0NkTkRkfD+(kR)ns<kT)%(kT) 
(kT-kR)n+(kR)R*  VT)D+(kT) 


D* \  dR 
R  kT  "  dn 


K  =  kr 


cp(x,z)  =  ^[A(H)e-nd<*>I  +  i“dK, 


B, 


t(x,z)  =  JrjB(x)e-n3Wz  +  i,lxdH,  x>0, 
B, 


(3.  33) 


(3.34) 


(3.35) 


and  where  denotes  the  branch-line  contour  in  the  upper  half- plane 
(see  Fig.  4).  The  full  field  is  thus  seen  to  have  the  required  form 
(2.  8)  for  x>0,  with  Afc  as  the  amplitude  of  the  transmitted  wave,  pro¬ 
vided  cp  and  iji  satisfy  the  radiation  conditions  (2.  5).  That  they  do  may 

f~ 2  2? 

be  proved  by  analyzing  the  asymptotic  behavior  for  large  Vx  +z  of 
the  integrals  in  (3.  35).  When  z  =  0,  such  an  analysis  may  be  carried 
out  with  the  aid  of  Watson's  lemma,  and  when  z/  0,  by  the  saddle 
point  method. 

When  x<0,  a  similar  procedure  applied  in  the  lower  half  of  the 
H-plane  leads  to  the  result 

‘Mx,  r.)  =  AQcpR+(x,z)  +  ArcpR_(x,  z)  +  cp(x,  z)  ,  | 


Y(x,z)  =  AQi|rR+(x,  z)  +  Ari(iR_(x,  z)  +  i|i(x,  z)  ,  x<0  , 


(3.  36) 
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where 


AQkR^ks"kd^D  ^kR^ns^"kR^ 
n+(kR)D‘(-kR)R'(-kR) 


(3.  37) 


R'<‘kR>  =  f 


H  =  -k 


R 


cp(x,z)  =  2“  jA(x)e“nd^H^z  +  1KxdH, 


B 


i|)(x,  7.)  =  J B(x)e  n 


s(*)z  +  iwd)(j  x<0  >, 


(3.  38) 


B. 


and  where  B_  denotes  the  branch-line  contour  in  the  lower  half-plane. 
Thus  the  full  field  has  the  required  form  for  x<0,  with  Ar  as  the 
amplitude  of  the  reflected  wave. 

Note  that  the  functions  cp(x,  z),  i)/(x,  z)  as  defined  by  (3.  35)  and 
(3.  38)  are  discontinuous  across  x  =  0,  z>0,  but  the  full  field  £  is 
twice-continuously  differentiable  in  the  interior  of  the  half- space. 

Also  note  that  had  (3.  32)  been  derived  formally  without 
recourse  to  the  assumption  of  complex  frequency,  the  inversion  con¬ 
tour  would  appear  to  pass  directly  through  the  poles  and  branch  points 
(since  they  would  now  lie  on  the  real  H-axis),  and  it  would  be  neces¬ 
sary  to  appeal  to  some  other  argument  to  determine  whether  this 
contour  should  be  considered  as  passing  above,  below,  or  through 
these  singularities.  The  assumption  of  complex  frequency,  in  addi¬ 
tion  to  providing  the  overlapping  domains  of  analyticity  required  for 
the  Wiener-Hopf  analysis,  serves  in  effect  to  resolve  this  dilemma. 
Having  served  its  purpose,  it  may  now  be  discarded. 


-36- 


It  is  now  desired  to  obtain  numerical  values  for  the  reflection 
and  transmission  coefficients  A^/A^,  A^/Aq.  All  quantities  in  (3.  34) 
and  (3.  37)  can  be  evaluated  straightforwardly  except  for  the  functions 
D+  and  D  .  It  is  proved  in  Appendix  B  that 

D  (-h)  =  — 7~ —  ,  Im{x)  >  -  6  ,  (3.39) 

D+(h) 

so  for  the  purpose  of  determining  A^Aq  and  A^/Aq  from  (3.34)  and 
(3.37),  it  is  only  necessary  to  compute  D^(k^)  and  D+(k^,).  This  is 
done  by  a  numerical  integration,  the  details  of  which  are  explained  in 
Appendix  B.  After  rewriting  (3.  34)  and  (3.  37)  in  dimensionless  form 
it  is  seen  that  the  coefficients  A^/Aq,  A^/Aq  are  functions  only  of 
the  dimensionless  layer  parameter  Nkg  and  Poisson's  ratio  v.  The 
dependence  on  Nkg  of  the  transmission  and  reflection  coefficients  is 
shown  for  v=l/4  in  Fig.  5,  The  coefficients  A^/Aq  and  A^/Aq  are 
complex,  but  for  the  range  of  Nkg  shown  the  imaginary  parts  never 
exceed  0.  4%  of  the  corresponding  real  part.  A  complex  coefficient 
simply  signifies  a  wave  whose  phase  is  shifted  with  respect  to  the 
incident  wave. 

A  check  of  these  numerical  results  may  be  obtained  by  com¬ 
puting,  from  (3.34)  and  (3.37),  the  asymptotic  forms  of  A^/Aq  and 
A j. / Aq  as  N-*0.  The  wave  number  k^,  is  a  function  of  N  defined 
implicitly  by  (1.  13): 

R[kT(N)]-Nk\^,(N)ns[kT(N)]  =  0  . 

This  equation  may  be  differentiated  with  respect  to  N,  yielding 
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R/[kT(N)]k,T(N)-k\^.(N)ns[kT(N)]-N^{k\5,(N)n  [k  (N)]}  =  0 


dN  s  T'  7  sL  T 


Upon  setting  N  =  0,  and  noting  from  (1.  13)  and  (1.  12)  that  k  (0)  =k 

T  R  * 

there  results 

kfk^n  (k„ ) 
k'  /oi  -•  s  R  R; 
kT'0)  “  R'(kR)  • 

NkVn  (k  ) 

kT  =  kR  +  —  f  *0(N1  as  N-0  . 

R 


Thus 


(3.  40) 


It  is  proved  in  Appendix  B  that  for  k>k 


,  /Nkf(k  +k) 

D  (k)  =j - =■-■%  -  +  o(Jjp)  as  N-0 

2(k  -k  ,) 
s  d 


Application  of  (3.  39)- (3.  41)  to  (3.34)  and  (3.37)  yields 
At  =  A0+  » 


(3.41) 


(3.  42a) 


Nk  k  n  (k  ) 

Ar  =  2R'(-k  )  +  as  N"° 

K 


(3.  42b) 


The  first  terms  of  (3.  42a)  and  (3.  42b)  are  plotted  in  Fig.  5  and  it  is 
clear  that  they  do  agree  with  the  numerical  computation  with  increas¬ 
ing  accuracy  as  N-0.  Equation  (3.  42a)  expresses  nothing  more  than 
the  fact  that  the  amplitude  of  the  transmitted  wave  approaches  that  of 
the  incident  wave  in  the  limit  of  vanishing  layer  thickness. 
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4.  Approximate  Solution  for  a  Very  Thin  Layer 

It  is  now  assumed  that  the  layer  is  thin  enough,  when  com¬ 
pared  with  the  wavelength  of  any  of  the  characteristic  surface  or  body 
waves,  that  the  layer  parameters  M  and  N,  defined  by  (1.  8),  may  be 
taken  as  perturbation  parameters  in  a  sense  to  be  described  pres¬ 
ently. 

Let  a  constant  a  be  defined  by 


a 


(4.  1) 


so  that 


M  =  aN  . 


(4.  2) 


The  perturbation  process  consists  of  allowing  the  layer  thickness 
2  h'  to  grow  from  zero  to  some  small,  non-zero  value  while  holding 
o  fixed.  Thus,  in  view  of  (1.  8),  (4.  1),  and  (4.  2),  the  problem  may 
be  expressed  in  terms  of  the  single  perturbation  parameter  N. 

The  differential  equation  (2.  6)  for  z)  remains.  By  using 
(4.  2),  the  boundary  conditions  (2.  7)  may  be  rewritten  as 


R£  =  0  ,  z  =  0  ,  -oo<x<  0  , 

[R  +  N(S  +  oT)ll  =  0,  z  =  0,  0<x<oo  . 


} 


(4.  3) 


The  assumed  form  (2.  8)  also  remains  appropriate,  with  the  wave 
number  k.^,  a  root  of  (1.  9),  but  now  it  is  assumed  that  the  amplitudes 
Ar,  Afc  appearing  in  (2.  8)  may  be  represented  by 


A  =  A(0)+NA(1)+  o(N),  > 

r  r  r 

=  A<0)+  N  a[1  ’  +  o(N)  asN-oJ 


(4.4) 
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and  that  the  radiated  wave  9  may  be  expressed  as 

?(x,a)  =  ^(0)(xf  z)  +  N2(1)(x,z)  +  o(N)  as  N-0  .  (4.  5) 

Each  term  on  the  right  of  (4.  5)  is  assumed  separately  to  satisfy  the 
radiation  condition  (2.  5). 

The  remainder  of  this  section  is  concerned  with  the  determi¬ 
nation  of  \  m  =  0,  1.  The  first  step  is  to  derive  the  equa¬ 
tions  governing  and  By  substitution  of  (2.  8)  into  (2.  6)  and 

{4.  3),  there  results  the  following  problem  for  <£: 


_L^  =^,  z>  0 ,  -oo<x<oo  ,  x  i  0, 

Rcp  =  0,  z  =  0,  -oo<x<0, 

{R  +  N(S  +  aJT)}  C£  =  0  ,  z  =  0  ,  0<x<oo  . 


J 


(4.6) 


The  requirement  that  _£  be  twice-continuously  c  ifferentiable,  when 
applied  to  (2.  8),  leads  to  the  jump  conditions 


W  -  Ao^R+(0-  z)  -  AA-(0'  z)  +  At^T+(0'  z)  =  0*^ 

[&x]  -  ikRAc^R+<0'  z)  +  ikRAr2R.(°.  a) 

+  ikTAtTT+(0,  z)  =  j0,  0<z<oo  , 

where 

fcp]  =  cp(0+,z)-cp(0",z)  . 


When  (4.  2)  is  substituted  into  (1.  9),  there  results  an  equation 
which  defines  kT  implicitly  as  a  function  of  N.  By  differentiating 
this  equation  with  respect  to  N  and  setting  N  =  0,  it  follows  that 


.  (1L  ^2 

T  r  dN 


N~0 


ir^)fkRkJna(kR)-a,kJ[nB{kR»+nd(fcR)]^  ■  (4‘  8) 
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Equation  (4.  8),  and  the  fact  that  =  kR  when  N  =  0,  lead  to 


kT  =  kR  +  Nk^  +  oW  as  N-0  . 


(4.  9) 


The  quantity  ^-p+(0,  z),  defined  by  (2.  8),  may  similarly  be 
expanded  in  a  power  series  in  N,  resulting  in 

cpT+(0,z)  =cpR+(0,  z)+NC£^(0,z)  +  o(N)  as  N-0,  (4.  10) 


where 


£^<o,  zh 


l-n<d1,ze-nd(kR>* 


|(P^)-PRns(1)z)e'ns(kR|z 


dN 


n(I)  .. 


N=0 


dN 


N=0 


Now  (4.  4),  {4.  5),  (4.  9),  and  (4.  10)  are  substituted  into  (4.  6) 
and  (4.  7),  and  the  coefficients  of  the  zeroth  and  first  powers  of  N  in 
each  equation  are  separately  set  equal  to  zero,  to  yield 

L  =  0  ,  z>0,  -oo  <x  <oo  ,  x^O; 

=  0,  z  =  0,  -oo<x<oo,  x£0; 

[cp(0)]-A09R+(0,z)-Aj>0)cPRJO,z)  +  Aj0)9R+(0>z)  =  0,  V,  (4.  n) 

^  -  ikRAo£R+<°-  zl  +  ikRAr0l2R-(°.  *> 

+  ikRAi0)2R+(°'z)  =  0f  z>  0  ; 


J 


-42- 


Ij  cp(  ^ )  =  o  ,  z  >  0  ,  - oo <x <cx) ,  x  i  0  ; 


j^cpn)_Q^  z  =  o,  - oo <x <0  ; 


-\ 


RfP(1>  =  -(S  +  a-T)C£(0) ,  z  =  0  ,  0<x<oo  , 

[cp{1)]  -  A^^  JO.zJ  +  A^^^.zJ  +  A^cp^^O^z)  =  0  A  (4.  12) 


[^>]  +  i^A^^JO,  z)  +ikRA{1  VR+(0,  z) 


ikRA(t0y^(0,z)+J:^^0)2R+(Ofz)  =  0,  z> 0  . 


J 


In  order  to  solve  (4.  11),  it  is  convenient  to  define  an  auxil- 
.  *(0) 

iary  function  cp  by 

$(°)  =  y°>_  Ao2R+H(x)  +  A^0)cpR_H(-x)  +  A(t°yR+H(x)  .  (4.13) 


When  (4.  13)  is  substituted  into  (4.  11),  the  equations  governing  9 
are  found  to  be 


■n(0) 


L  <£'  =  () ,  z> 0  ,  -oo <x<oo  , 

z  =  0  ,  -  oo  <x  <oo  , 
[$(0)]  =  [cp(,0)]  =  0,  z>0  .  J 


(4.  14) 


■—  x 


\(0) 


In  particular,  note  that  the  differential  equation  for  ^Pv  '  holds 

throughout  the  interior  of  the  half- space.  From  (4.  13)  it  is  seen  that 

-  (0) 

cp'  is  everywhere  outwardly  propagating.  This  fact  together  with 
(4.  14)  implies  that  Thus,  from  (4.  13), 


cp(0)  =  <A0-  A[0))cpR+H{x)  -  aJ.0)C£r_H(-x) 


(4.15) 


But  since  cp^  is  a  radiated  wave,  it  cannot  contain  any  surface 


waves,  so  from  (4.  15), 
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A 


(°)_  A  A(0)- 
t  “  A0  1  r  " 

C£(0)(x,  z)bO  . 


(4.  16) 


Equations  {4.  16)  simply  express  the  fact  that  in  the  limit  as  the 
layer  vanishes,  the  full  field  consists  only  of  a  positive-traveling 
Rayleigh  wave  of  amplitude  Aq. 

In  view  of  (4.  16),  Eqs.  (4.  12)  become 


L  <P^  =  £,  z>0,  -oo  <x<oo  ,  x  /  0  ; 

R  ^  =  0_ ,  z=0,  -oo<x<oo  ,  x/0; 
fcp(1)]  .  aJ.1)^r_(0,z)  +  a[1)C£r+(0,x)+A0^(^(0,z)  =  0 
+  ikRA(rJ  ^.(0,  z)  +  ikRA{1)£R+(0,  z) 

+  ikRA(t°  V^(0,  z)+ik(^)A(t0)cPR+(0,  z)  =  0  ,  z>  0 


(4.  17) 


The  solution  of  (4.  17)  is  facilitated  by  the  introduction  of  an 
"'(1) 

auxiliary  function  defined  by 

^(1)=  cp(1)+ A^'cpR  H(-x)  +  A(tJ)^R+H(x)  +  A^^O,  z)eikRxH(x) 


A  k;., 

+  [^R+H(x)  +  cpR+(0,z)e“lkRXH{-x)J  .  (4.  18) 

A  (  1  ) 

When  (4.  18)  is  substituted  into  (4.  17),  the  equations  for  cp  are 
found  to  be 

L  ^ 1  ^  =  AqH(x)L  z)e^kRx]  ,  z>0,  -oo<x<oo,  (4.19a) 


A  k(i) 

^ifeR+(°.-)e-ikRX}H(-X) 

+  A0R{cp(^(0,z)eikRx}H(x),  z=0, 


-00<X<00  , 


(4.  19b) 
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[£(1>]  =  [i!^]  =  o,  z>0  .  (4.19c) 

a/Q)  .  , 

As  in  the  case  of  cp'  ,  the  function  0'  '  is  smooth  enough  that  the 
differential  equation  (4.  19a)  holds  throughout  the  interior  of  the  half¬ 
space. 

The  right-hand  sides  of  (4.  19a,  b)  may  be  simplified  some¬ 
what.  Note  that 

cpT+(x,z)  =  ^T+(0,z)elkTX 

=  f"c£R+(0,  2)  +  Nqp(J|(0f  z)  +  0(N2)]  elkRX  fl  +  iNk<? }x  +  OfN2)"1 
=  9>R++N^fr1](0,z)e1KRx+ik^)xqpR+]+0(N2)asN-*0  .  (4.20) 

Since  h  J£<j,+(x,  z)  ^  0,  each  term  of  the  latter  form  of  (4.  20),  when 
operated  upon  by  L,  must  vanish.  Thus 

L{^|{0,z)eikRX}=  -ik(^JJx2R+(x,z)}  .  (4.21) 

From  (1.  8),  (1.  2),  (1. 4),  (2.  7),  and  (4.  2)  it  follows  that 

{R  +  N(S  +  aT)}cpT+=  0,  z  =  0  .  (4.22) 

Substitution  of  (4.  20)  into  (4.  22)  leads  to 

R^R+  +  4i^T+{0'z)elkRX^  +ikT  =fx^R+^  +  (I  +  °1)-R+] 

+  0(N2)  =  0  as  N-0  ,  z  =  0  .  (4.  23) 

Each  term  of  (4.  23)  must  likewise  vanish  separately  so 

RtjP^fO,  z)eikRX1  =  -ik!r1)R{xCPR+3  -  (S  +  aT)fR+  ,  z  =0  .  (4.  24) 

The  right-hand  sides  of  (4.  21)  and  (4.  24)  may  be  easily  computed  by 
using  Leibnitz'  rule,  and  by  noting  that  L  ^PR^.  =  0  and  that 
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H^r+  =  £  at  z  =  0.  The  first  term  in  the  right-hand  side  of  (4.  19b) 
may  also  be  computed  directly.  These  results  are  then  incorporated 
into  (4.  19)  to  yield 


£$<1)=2A0kRl41,;eR+H(x),  Z>a,  -oo <x<co 


(4.  25a) 


RCp'  '= 
where 


,d)_  A0kT><24-^,  ...  .|r 

- 2r-h(-x)  +  V 


ll  ik_x 


e  R  ,  z  =  0,  -oo<x<oo, (4.  25b) 


al  ' 


kT1,<24-kS2) 


2  2 
2k£-k 
R  s 


a2-i[^  +  WU,.g]. 


R 

2 


nd(kR}"2n  (k  S)]» 

S  xv 

2 


(4.26) 


J 


Equations  (4.  25)  may  be  solved  with  the  aid  of  the  Fourier 
transform  (3.  7).  It  is  once  again  assumed  that  the  frequency  ou  has  a 
small,  positive,  imaginary  part.  The  transform  of  (4.  25a)  is  then 
found  to  be  ^ 

«) .  a^i>(K. .) , 


»2i(D 


-n^nfflVz)  a 

dz *  s  »(*- 


i(x-kR) 
2Ankck(^)(3pe*ns(kR)z 


>  (4.27) 


_ R_ 

i(x-kR) 


z>0  ,  -6<Im{x}<6  ,  J 

where  6  is  a  constant  such  that  0<6<Im{kd}.  The  solution  of  the 
homogeneous  version  of  (4.  27)  is  given  in  (3.  11),  and  again  A '(h), 
B'(k)  must  vanish.  The  general  solution  of  (4.  27)  is  thus 
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2A„k„k,i)Vnd,kR>* 


i(H-kR)  (n+kR) 


2A„k„k,’lp_o'”slkR,z>  (4.28! 

*(1)(k,Z)=B(K)e-n3Wz-— - ,  | 


i(H-kR)  (x+kR) 


z>  0  ,  -6<Xm{K]  <  6  . 


J 


Application  of  the  transform  (3.  7)  to  the  boundary  conditions 
(4.  25b),  followed  by  use  of  (4.  28),  yields 


/2H2-k2 

s 


2inn  (x)\  A(h) 


s 

2  ,2 


2AOkRkT* 


!kR-k2 

R  s 


-2ixnsncR)) 


2ixnd(K)  -(2k  -ksj/  B(„)  i(H-kR)^(H4kR)  \-2iKnd(kR)  -(2kR-k^)/  k 


AQk^)(2k2-k2)  r  1 


R  s' 


ikR(n+kR) 


<  Im  {*]  <5  .  (4.29) 


This  is  a  system  of  two  simultaneous,  linear  equations  in  the  two 
unknowns  A(h),  B(h).  It  possesses  a  unique  solution  provided  the 
determinant  of  the  coefficient  matrix  does  not  vanish.  That  determi¬ 
nant  is  found  to  be  -R(h),  where  R(h)  is  given  by  (3.  17b)  and  is 
known  to  vanish  only  at  h  -  ±  kR.  Since  these  points  lie  outside  the 
domain  of  validity  of  (4.  29),  A(h)  and  B(x)  may  be  computed  imme¬ 
diately,  and  are 
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2ArtkJc' 


,  .  (2n2-k2)(2k2-k2)  ^ 

A<* : 1  *  ■ -  .o-^TO  { - % - +  4“S'’‘>”d<kR>} 


+  T 


AOk<T(24-k'> 


ikR(n+kR)R(H) 


r  9  ?  Hn  (h)  ■, 

{-(2»  -k')-(2kR-ks,k^)} 


0 


i(H-kR)R{K) 


|-aj(2K2-k2)  +  2mng(K)a2|  , 


(4.  30a) 


«<»> = i(,-X°K£ro  K<  v<2-k'» + 2i“d<kR)4'24-kf »} 


.(i) 


.rM-'-R 


A-kfJ,\2k2-k2)  ,  2ik  n  (k  )  ,  2 

+  .,0  J  I,  (2iKn  .(H)  +  fi-l?A.(2K2-k2)} 

ikR(H+kR)R(K)l  d  2k;-k2  sJ 

R  s 


i(K-kR)RU) 


|2iHn^(K)a ^  +  a2{2H2-k2)|  ,  -6<Im{x}<&  .  (4.30b) 


Equation  (4.  30a),  for  example,  may  now  be  used  in  the  inversion 
formula 

OO 

$(1)(x,z)  =  2^  J[A(*)e"nd 


2A  k  kliVnd^R^z-,  . 

<»>»-  O  R. 2l - ]e,Rxd».  (4.311 


-CD 


i(n-kR)  (x+kR) 


which  follows  from  (3.  8)  and  (4.  28). 

The  singularities  of  the  integrand  in  (4.31)  are  branch  points 

at  ±k  ±k  .  and  poles  at  ±k_,.  The  contributions  due  to  the  residues 
d  s  K 

at  the  poles  represent,  as  in  the  previous  section,  the  surface  waves. 
For  x>  0  the  contour  may  be  closed  in  the  upper  half-plane,  as  shown 
in  Fig.  4,  so  in  this  case  it  is  the  pole  at  kR  which  is  of  interest. 
Four  separate  terms  appear  —  three  due  to  A(k)  and  the  fourth 
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arising  from  the  last  term  in  (4.  31).  The  residue  of  each  term  may 
be  computed  separately.  The  first,  third  and  fourth  terms  have 
double  poles;  the  second,  a  simple  pole. 

The  double  poles  will  lead  to  secular  terms,  i.  e.  ,  terms  con¬ 
taining  x  or  z  as  a  factor.  It  will  be  seen,  however,  that  none  of 
these  terms  appears  in  the  asymptotic  expansion  of  *P^(x,  z)  as 
x-*±oo.  This  is  necessary  in  order  that  the  asymptotic  expansion  be 
valid;  otherwise  the  representation  of  V^(x,  z)  would  not  be  uni¬ 
formly  valid  in  x  and  the  asymptotic  expansion  could  not  be  justified. 
In  fact,  it  is  this  very  requirement  which  necessitated  the  seemingly 
roundabout  formulation  employed  in  this  section. 

From  (4.  18),  for  x>0, 


cp^(x.z)  =  cp^(x,  z)  -  A^cpR+(x,  z) 

A  k(1) 

(  I  )  ilr  v  0  T 

-  A0<p^(0,z)e1KRX-  <pR+(x,  z) ,  x>0  .  (4.32) 

When  the  four  residue  contributions  from  are  combined  with  the 

latter  three  terms  of  (4.32),  there  results 

cp^^(x,  z)  ~  (A^  -  A^)  qpR+(x,  z)  as  x-»oo  ,  (4.33) 


where 


AnR"(kp)k  f 

A  =  - f  rkRns(kR^  +  ®  ^ns^kR)  +  nd^kR^J 

1  2[R'(kR)]  1  R  s  R  s  R  d  R  J 


A0ka2(2ki-k?) 


ng(kR)R  ,~R 


R  sT.  Za 

rWniR'  kD  J* 


(4.  34) 
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Note  that  the  secular  terms  due  to  the  double  poles  have  been  can¬ 
celled  by  each  other  and  by  the  third  term  on  the  right-hand  side  of 
(4.32).  Now  since  <P^(x,  z)  must  obey  the  radiation  condition  (2.  5), 
it  cannot  contain  any  surface  waves,  so  from  (4.  33)  it  follows  that 


(4.  3  5) 


The  coefficient  may  be  determined  by  a  procedure  analo¬ 

gous  to  that  described  above,  but  with  the  integration  contour  closed 
in  l°wer  half  K-plane.  Because  A(h)  has  only  a  simple  pole  at 

=  the  computation  is  somewhat  simpler.  The  final  result  is 


A 

Ar,=  2k“R'(-kR){-kRns,kR)ks  +  i,'ks1;ns(kR)-nd<kR)]}-  <4'  36> 


When  (4.  16)  and  (4.  34)  -  (4.  36)  are  used  in  (4.  4),  the  trans¬ 
mission  and  reflection  coefficients  may  be  computed  to  first  order  in 
N.  The  resulting  approximations  to  the  coefficients  are  shown  in 
Fig.  6  for  v  =  1/4. 

If  the  layer  is  massless,  M  =  0,  a  =  0,  and  the  results  may  be 
compared  with  those  of  the  previous  section.  This  is  done  in  Fig.  5 
for  v  =  1  /4,  and  the  re  is  seen  to  be  good  agreement  as  N  -*  0. 


a 

o 


Figure  6.  Approximate  Transmission  and  Reflection 
Coefficients  for  a  Very  Thin  Layer. 
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5.  Conclusions  and  Comments 

The  problem  of  the  scattering  of  a  Rayleigh  wave  by  the  edge 
of  a  thin  surface  layer  is  amenable  to  solution  by  the  Wiener-Hopf 
technique  in  the  case  where  che  density  of  the  layer  vanishes.  The 
amplitudes  of  the  reflected  and  transmitted  surface  waves  computed 
by  such  an  analysis  are  shown  in  Fig.  5  for  v  =  1/4. 

A  perturbation  process  leads  to  an  approximate  solution  for 
the  case  of  a  layer  possessing  both  mass  and  stiffness.  This  analysis 
yields  the  transmission  and  reflection  coefficients  shown  in  Fig.  6, 
for  v=  1/4.  The  results  of  the  two  analyses  agree  closely  for  very 
thin,  massless  layers. 

It  is  reasonable  to  ask  whether  the  results  of  Section  3  might 
be  extended  to  the  case  of  a  layer  having  mass,  by  treating  the 
layer's  mass  parameter  M  as  a  perturbation  parameter,  while 
allowing  the  stiffness  parameter  N  to  remain  arbitrary.  Such  a  pro¬ 
cess  leads  to  a  sequence  of  problems,  the  zeroth  order  problem 
being  precisely  that  solved  in  Section  3.  This  procedure  has  been 
attempted,  but  the  computations  required  to  obtain  the  first  order 
solution  were  found  to  be  of  prohibitive  complexity.  Part  of  the  diffi¬ 
culty  arises  from  the  fact  that  the  zeroth-order  solution,  which  is 
known  only  as  a  Fourier  inversion  integral,  appears  as  a  forcing 
function  in  the  first-order  problem.  This  necessitates  the  factoriza¬ 
tion,  by  the  formal  procedure  used  in  Section  3,  of  quantities  contain¬ 
ing  the  zeroth  order  solution. 
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Appendix  A  -  Roots  of  R  (a) 

In  this  section  the  roots  in  a  suitably  cut  complex  plane  of 
the  function 

R*(k)  =  (ZHZ-kZ)Z  -  4k  2n  (x)n.(x)-NkVn  (h),  N>0,  (Al) 

'  s  s  Q  s  s 

are  located. 

It  is  assumed  that  uu  =  0,  so  that  kg  and  kd  are  real.  The 
branch  cuts  in  the  complex  n-plane  are  taken  as  shown  in  Fig.  7, 
except  that  the  radii  of  the  cuts  in  the  vicinity  of  the  origin  is  vanish¬ 
ingly  small.  The  functions  ns(H)  and  n^H.)  are  then  defined  by 
requiring  that  they  be  positive  for  H>kg.  When  ng(?t)  an<3  nd(n)  are 

continued  analytically  throughout  the  cut  H-plane,  it  can  be  seen  that 

•>:  $  >!« 

R  00  is  even  in  n  ,  i.  e.  ,  R  (h)  =  R  (-k).  Therefore  it  is  only  neces¬ 
sary  to  consider  the  right  half-plane  Re{n}a0. 

The  number  of  roots  in  the  right  half-plane  may  be  determined 
by  the  principle  of  the  argument  (see  pp.  271-272  of  [22]).  The 
principle  may  be  stated  as  follows:  Let  a  function  be  analytic  inside 
and  on  a  closed  contour  C  except  for  at  most  a  finite  number  of 
poles  inside  C.  Let  f  have  no  zeros  on  C  and  at  most  a  finite 
number  of  zeros  inside  C.  Then 

No- V^Ac[argfW]  (A2> 

where  Nq  is  the  number  of  zeros  of  f  inside  C,  a  zero  of  order 
mQ  being  counted  mQ  times,  is  the  number  of  poles  of  f  inside 
C,  a  pole  of  order  being  counted  times,  and  A^[argf(«.)]  is 
the  change  in  the  argument  of  f (h)  as  C  is  described  in  a  positive 


sense. 
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A1  igure  7.  Complex  H-Plane  for  Determination 
of  Zeros  of  R*(h). 
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The  curve  C  is  taken  as  shown  in  Fig.  7.  It  is  clear  by  in- 

spection  of  (Al)  that  R  (x)  has  no  poles  inside  or  on  C,  so  (A2), 

>!<  £ 

when  applied  to  R  (n),  will  give  the  number  of  zeros  of  R^(h)  inside 

>|(  }IC 

C,  provided  that  R  (x)  has  no  zeros  on  C.  To  pro^e  that  R  (h)  ^0 
on  C,  each  segment  of  C  must  be  examined  separately.  The  appro¬ 
priate  values  of  ng(x)  and  nd(x)  for  each  segment  of  C  are  deter¬ 
mined  by  analytic  continuation  from  the  half-line  Im{x}  =0  , 
Re{xl>kg.  These  and  the  corresponding  formulas  for  R*(x)  are 
given  in  Table  1  in  terms  of  5  and  r|,  where  h  =  §  +  it).  By  referring 
to  the  table,  the  following  facts  are  noted:  On  C 2  and  C^, 

Im{R  (x)}  =0  only  at  ti  =  0,  where  Re {R  (h)]  /  0.  On  C 3  and  , 

hn {R  (h)  }  =  0  only  at  5  =  0,  where  RefR'  (x)]  t  0  On  C.  and  Cc, 

4  D 

Re{R(x)}=0  only  at  §  =  k  /«/ 2~,  where  Im{R^(x)}  /  0.  Hence 

s 

r"1  (x)  / 0  on  C2  -  C?. 

i  0  » l 

On  Cj,  ng(x)~nd(x)~re  as  r-*oo,  where  8  =  Tan”  (r|/5),  so 
from  (Al), 


R*(x)~[2(re10)2-  k2]2- 4(rei0)2(rei0)2- Nk2(rei0)2reifl 
s  s 

~  ~  Nk2r^e0*0  as  r-*oo  (A3) 

s 

3fC 

Thus  R  (h)  ^0  on  C,  so  the  principle  of  the  argument  may  be  applied. 

One  way  to  determine  (l/?/n-)Ac[argR  (x)]  is  to  plot  the 

.  #  # 

trajectory  of  R  (x)  in  the  complex  R  -plane  as  x  describes  C  in  the 
complex  x-plane,  and  then  simply  count  the  net  number  of  times  the 
origin  is  encircled  in  a  counterclockwise  direction.  The  complex 

V  _ 

R  -plane  is  shown  in  Fig.  8.  Point  A,  which  corresponds  to  x  =  -ir, 
is  taken  as  a  starting  point.  Its  image  A'  in  the  R  -plane,  by  either 


Segment  Range  of 
of  C 


Tabic  1  -  Formulas  for  n  (h),  n  ,(k),  and  R  (a). 


f . 

1 

$ 

i 

L 

i 

i 

% 
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II 

^  2  3 

(A3)  or  Table  Al,  corresponds  to  R  rw-iNkgr  .  As  is  de¬ 
scribed,  0  increases  by  it  while  r  remains  fixed,  so  by  (A3), 
arg(R*)  increases  by  3tt  while  |R  |  remains  constant  and  large.  As 

c  is  described,  Re{R^}  never  vanishes  since  here  Re{R  (k)}  =R(h), 

2 

and  R(h)  is  known  not  to  vanish  on  the  imaginary  axis  (see  p.  190  of 

[16]).  Thus  on  Re{R*}>0  and  ImfR'^O,  equality  holding  only 

at  r|  =  0.  On  C3>  Re{R>0  and  Im{R  }  50,  equality  holding  only  at 

P  =  0.  On  C  ,  Re  [R*]  =  0  only  at  %  =  k  If? ,  where  Im{R}<0. 

4 

Otherwise  Re{R*l>0  and  Im{R!]  £0,  equality  holding  only  at  §=kg. 
The  images  in  the  R*  - plane  of  the  remaining  segments  C5,  C&,  and 

C?  are  the  complex  conjugates  of  those  of  C^,  and  C^,  respec- 

U  * 

tively.  Upon  returning  to  the  starting  point  A'  in  the  R  -plane,  the 

| 

origin  has  been  encircled  once  in  a  counterclockwise  direction,  so 
by  the  principle  of  the  argument,  there  must  be  precisely  one  root 
of  R" (k)  in  the  right  half  of  the  K-plane,  and  it  must  be  simple. 

To  locate  this  root,  consider  R  (k)  for  K  =  ?,  kg=5§£kR.  ^ote 
from  (A  1)  that  R*(kg)  =  kg  >  0,  and  R* (kR)  =  -  NkgkVkR-kg  • <  0 
(since  R*(r)  =  R(h)  -  NkgR2ng(x)  and  R(kR)=0).  Hence  the  root,  which 
is  denoted  by  kT,  must  lie  in  the  interval  (hg,kR).  Since  R.  (h)  is 
even  in  h,  it  has  another  root  at  -kT> 

1 

! 

1 
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Appendix  B  -  Factorization  of  D(h);  Properties  and  Evaluation 
of  D+(*)  and  D"(h). 

The  function  D(k)  is  defined  as 
2(k2-k2)R*(H) 

D(h)  =  - 1 - - -  f  -6*Im[K}s6f  (Bl) 

Nk^n  (h)R(k) 

b  b 

where 

R(«)  =  (2m2-  k2)2  -  4«.2ng(H)nd(K)  , 

R*(h)  =  R(k)  -Nk*K2n  (k), 

!~?  71  r~2  T 

n  (k)  =  Vh  -  k  ,  n  (k)  =  Vk  -  k  ,  , 
s  s  u  a 

and  where  6  is  a  constant  such  that  0<6<Im(kd}.  The  functions 
ng(n)  and  nd(n)  are  rendered  single- valued  by  taking  branch  cuts  as 
shown  in  Fig.  4  and  by  requiring  both  to  approach  too  as  M.-*oo  along 
the  positive  real  n-axis. 

& 

The  functions  R(h)  and  R  (h)  have  the  following  asymptotic 


fo  rms : 
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R(h)  =  4k 


4k 


2  2 


2  1 


k-3H 


■2K2(lvkd)  +  0(1) 


■2|H|2(k^-kp  +  O(H) 


TT1 

2  o  /  k 

R  '"(h)  =  R(K)-NkgH7  1--| 

K 


Nk2H3+  0{  |k|) 
s 


W) 


=  -Nk2  U|3+  0(  |h.|*0  as  |h|-*ot,  -6£lm{x}^f),  Re{H}>O.J 


In  (B2)  use  was  made  of  the  fact  k=|h|+0{1)  as  |x|-*oo,  -6^Imfn}^6, 
RefH}>0.  Since  R(h)  and  R  (h)  are  even  functions  of  k,  the  esti¬ 
mates  (B2)  also  hold  as  |h|  -*oo,  -6^Im{H}^6,  Re[x,}  <0.  Now  if  (B2) 


is  substituted  into  (Bl)  it  is  seen  that 

D(h)  =  i  +  0(  H"1)  as  |h|^oo  ,  -6^imfK]^6  . 


(B3) 


The  functions  D+(k)  and  D~(h)  are  formally  defined  as 

oo^i^ 

D  (h)  =  exp  {"  2^1  J  -5<IxnM<t  ,  (B4) 


-OOTiS 


where,  in  view  of  (B3),  LogD(a)  is  rendered  single- valued  by  requir¬ 
ing  it  to  approach  0  as  jty| -Hoo ,  Re{o-}>0.  The  function 

D(a)  has  branch  points  at  ±kg  and  ±k^,  simple  poles  at  ±k^,  and 
simple  zeros  at  ±kr^  (see  (Bl)}.  Since  none  of  these  points  lies 


t 


-60- 


within  the  strip  - S^trn f  a]  £6,  the  function  Log  D(a)  is  analytic 
throughout  the  strip.  Since  D(o)  is  an  even  function  of  a,  it  must 
also  be  true  that  LogD(a)-*0  as  |or|-*oo,  -5^Im{o}^&,  Re{o}<0.  In 
view  of  (B3) 

Log  D(a)  =  Odcvl  *)  as  |a[-*oo,  -  6£lm  [a]  ^6  ,  (B5) 

so  the  integral  appearing  in  (B4)  is  convergent. 

Now  by  using  (B4), 


co+i6  oo-i6 

=  eXp{--L  r  i2£Dkldff+  *  LogDiS-W} 

rL  2tu  J  a-K  2rriJ  ar-x  J 


-oo+i6 

.oo-ift  -oo+ifi, 


-oo-i6 


(B6) 


The  only  singularity  of  the  integrand  in  (B6)  between  the  two  infinite 
contours  is  a  simple  pole  at  a=K.  The  two  contours  may  be  joined 
by  two  vertical  segments  of  length  26  at  Re{»}  =±oo,  and  the  residue 
theorem  applied  to  the  resulting  closed  contour  to  yield 


=  exp { Log  D(x)3  =  D(h)  ,  - 6 <Im { x }  <  6  .  (B7) 

D+(H) 

The  "factorization"  procedure  exemplified  by  (B4),  (B7)  follows  the 
method  described  by  Noble  (p.  15  of  [21]). 

A  theorem  concerning  the  domain  of  analyticity  of  a  function 
defined  by  an  infinite  integral  (see  p.  110  of  [24]),  when  applied  to 
(B4),  leads  to  the  result  that  D  (x)  is  analytic  for  Im(x1>-6,  and 
D  (x)  is  analytic  for  Im{x]<6. 
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The  singularities  of  the  analytic  continuation  of  D  (h)  into 
the  lower  half-plane  are  determined  by  combining  (B7)  and  (Bl)  to 
yield 

Nkfn  (h)R(h) 

D+{k)  =  D“(H) - ■#■■■%■  > -  . 

2(k*-k*)R  (*) 


Thus  D  (h)  has  a  simple  pole  at  -kT  and  branch  points  at  -k  and 


■k^.  Conversely, 


+  2(k*-k?)R*(>0 

D"(h)  =  D  (K)  - ^ - - - 

Nk  n  (h)R(x) 

s  s 


so  the  singularities  of  the  analytic  continuation  of  D"(x)  into  the 
upper  half-plane  are  a  simple  pole  at  k^  and  branch  points  at  kg 
and  k 

CL 

+ 

It  is  clear  from  (B4)  that  D  (x)  and  D  (x)  can  have  no  zeros 

in  their  respective  domains  of  analyticity.  It  also  follows  from  (B4) 

and  the  evenness  of  D(cv)  that 

oo+i6  oo-i& 

exp{-  J 

-co+ifi  -oo-i6 

oo-i6 

=  expf^n  J  ~gS-^~'dQf}  =  D+(*>,  im{x}  >-6  ■ 

-oo-ifi 

Now  it  is  desired  to  show  that  D*(x)-*1  as  |k|-*oo, 

If  the  integral  appearing  in  (B4)  for  D*(x)  is  defined  as  1(h),  it  is 
sufficient  to  show  that  I(x)-«0  as  |x|-«oo,  Im{x}>-6.  In  view  of  (B5) 
and  the  continuity  of  Log  D(a),  there  exists  a  positive  constant  P 
such  that  |Log  D(a,)|<  P|o|  -6£lm  fa)  £6.  Therefore 
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where  0  =  arg(n).  It  can  be  shown  that  the  integral  appearing  in  (B8) 
is  0(Log|a|)  as  |x]-*oo  ,  Im{Hl  > -6,  so  I(k)  =  0(  |h.|"  1  Log  |  k|) -•  0  and 
D  (h)-*1  as  |h|  -too, 

The  quantities  D+(k„)  and  D^{k_)  must  be  evaluated  numeri- 
cally  in  order  to  compute  the  reflection  and  transmission  coefficients. 
The  integral  1(h)  appearing  in  (B4)  for  D+(h)  may  be  evaluated  by 
choosing  branch  cuts  as  shown  in  Fig.  9,  closing  the  contour  in  the 

lower  half-plane,  and  allowing  Im[uu]  to  vanish.  Note  that  the  zeros 

.1. 

'■N 

of  R(a)  and  R  (cv)  give  rise  to  logarithmic  branch  points  at 
“=±kR.  The  integrals  along  the  large  quarter-circular  arcs 

vanish  as  R-*oo,  leaving  only  the  integral  along  the  branch-line  con¬ 
tour.  This  contour  is  subdivided  into  ten  segments  numbered  sequen¬ 
tially  (see  Fig.  9).  The  integral  along  each  segment  may  be  combined 
conveniently  with  that  from  the  corresponding  segment  immediately 
across  the  branch  cut,  yielding 
oo-i6 

«x>  e  |  =  I1+10+  I2+9+  IJ+8+  I4+7+  I5+6.,B9) 

-oo-i6 


(The  contribution  from  the  small  semi-circles  about  the  singularities 
vanishes).  By  using  the  appropriate  complex  values  of  the  square 
roots  and  the  logarithm,  the  following  formulas  may  be  derived: 
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co 


1  + 


10=  2iJ 


Tan 


-1 


(2p2+l)2-4p2^ 

'ATv 

n-zi 

p  +Y 

dp 

^  Nkgp2  Vr)2+ 1 

- 

p-in/k  ’ 
s 

(BlOa) 


*2+9  "  2l 


Tan 


-1 


0 

1 


2  "z-  -2/i"-?"7y2  -Z 


(2g  -ir*4g  Vl-g^VY  -g 

Nk  gVl-g2' 

s 


d£ 


g+H/kg  ' 


3+8 


=  21J{r 


Tan 


-1 


wn  77  J 


RelE(g) 


*1 


g+K/k 


(BlOb) 


(BlOc) 


where 


E(g)  = 


-Nk^Vl-g2-  4g2\/l-^2Vg2-Y?-  i(2g2-  l)2 
(2g2-!)2-  4ie2/l-F2/?2- V2 


*4+7  "  0  ' 


J5+6  =  -2lTiL°g 


/kR+n\ 

Imv’ 


(BlOd) 

(BlOe) 


-1, 


where  Y^^/k^  and  where  -Tr/2sTan  (u)£tt/2,  -oo<u<oo.  The 

integrals  (BlOa  -  e)  have  been  evaluated  numerically  with  n  =  k„  or 

K  =kT>  then  the  results  used  in  (B9)  and  (F4)  to  give  the  numerical 

+  + 

values  of  D  (kR)  and  D  (k,p).  Note  that  since  y  depends  only  on 
Poisson's  ratio  v,  D^(kR)  and  D^(k^,)  depend  only  on  v  and  Nkg. 

The  asymptotic  form  of  D+(kt>)  or  D+(k„)  as  N-*0  is 
obtained  by  asymptotically  evaluating  (BlOa-e)  as  N  -  0.  The 
results  of  this  lengthy  analysis  are 


-65- 


5+10=  •  4  +0<NLogN), 

ry+>c/k  I 

X2+9°  -riLog^-^^  +  °(1|, 
u  s 

rl+K/k  -I 

I3+8=  -«L°g|_tli/FJ+  °(1), 

^5+6  =  ^(N)  as  N-*0,  kd<K<oo  . 

When  combined,  the  estimates  (Bll)  yield 


(Bll) 


2  2 

1(H)  =*i^g[fit^i£c7ig]-  ^T  +  odl  as  N-  0  . 


This  is  substituted  into  (B4)  to  yield 


/n(FTmF  .tt 

D  (k)  ~  / - - — =-  e  14  as  N-*0,  k  <K<oo  . 

V  i()-Y2) 
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Appendix  C  -  Certain  Aspects  of  the  Solution  in  the 
Case  of  a  Massless  Layer. 

It  will  now  be  proved  that  the  formal  solution  (3.  32)  satisfies 
the  differential  equation  (3.4)  for  -oo<x<oo,  z>0.  The  steps  will  be 

A 

carried  out  for  cp;  the  procedure  for  ij;  is  completely  analogous. 

Appeal  is  made  to  the  following  theorem  (see  p.  59  of  T24]): 
Suppose 

(a)  there  exists  an  e>0  such  that  f(x,ft)  and  8f/3x  are  continu¬ 
ous  for  x  -  e<x<x  +  e  ,  -Ks  k<  K,  for  all  K>0, 

o  o 

CO 

(b)  f(x,  K)dx  is  convergent  for  x  -  e<x<x  +  e,  and 

J  o  o 

-00 

oo 

r  3f 

(c)  I  -g^dH  is  uniformly  convergent  for  xQ-e<x<xo+e.  Then 
“00 

oo  oo 

Jf(x,*)d*  =  Jijd*  at  x  =  xo. 


-oo 


-oo 


Let  z  =zq>0  and  hold  it  fixed.  Apply  the  above  theorem  to 


the  function 


f(x,X»  =  AM<fndWVH,'X 


(Cl) 


By  (3.  31),  (a)  is  satisfied  for  any  e>0.  From  (3.  31)  and  (Cl), 

3  I  I 

f(x,  h)  =  0(  |h|  ^  e  !*'Zo)  uniformly  in  x,  and  f(x,  k)  is  bounded  for 
-oo<h<oo,  so  (b)  is  satisfied  To  verify  (c),  note  that 
df/3x  =  0(  |  h|  ^e~  ^z°)  uniformly  in  x  and  that  3f/0x  is  bounded  for 
-oo<n <oo,  so  there  exists  a  constant  P  such  that 


oo 


Bf_ 

9x 
-oo 


I 


d« 


oo 


fiP 


J*|?l|^e  ^  Z°dH<oo  ,  x_-e<x<x^+G  . 


( C  2 ) 


-oo 
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Thus  by  the  M-test  the  integral  is  uniformly  convergent  and  (c)  is 
satisfied.  Hence  the  theorem  holds,  and  the  differentiation  with 
respect  to  x  may  be  carried  inside  the  integral  for  -oo<x<oo,  z>  0. 

Due  to  the  exponential  factor  in  the  integrand,  the  above  process  may 
be  repeated  for  any  number  of  subsequent  derivatives  with  respect  to 
x  in  the  interior  of  the  half-space.  A  similar  procedure  may  be  ap¬ 
plied  for  derivatives  with  respect  to  z  and  for  cross-derivatives.  Thus 
derivatives  of  cp(x,  z)  of  all  orders  exist  in  the  interior  of  the  half¬ 
space  and  may  be  computed  by  formal  differention  inside  the  integrals. 

It  may  now  be  immediately  verified  that  the  differential 
equation  (3.  4)  is  satisfied  by  taking  derivatives  inside  the  integrals 
in  (3.  32). 

To  prove  that  the  derivatives  of  cp(x,  z)  are  continuous,  use 
is  made  of  the  following  theorem  (see  p.  54  of  [25]): 

If  a  function  of  g(x,  z)  has  uniformly  bounded  partial  deriva¬ 
tives  in  a  region  R,  then  g(x,  y)  is  continuous  in  R. 

To  apply  this  to  cp(x,  z),  consider  tp(x,  z)  in  any  closed, 
bounded  subdomain  of  the  half- space  -oo<x<oo,  z>0.  The  estimate 
(C2)  proves  that  9cp/8x  is  uniformly  bounded,  since  the  constant  P 
must  be  independent  of  zq  (see  (Cl)).  Correspondingly,  9<p/9z  is 
uniformly  bounded,  so  the  conditions  of  the  theorem  are  met  on  any 
closed,  bounded  subdomain  of  the  half- space. 

The  above  argument  proves  that  cp(x,  z)  is  continuous;  it  may 
be  repeated  to  prove  that  all  derivatives  of  $(x,  z)  are  likewise  con¬ 
tinuous  in  the  interior  of  the  half-space. 
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The  procedure  for  verifying  the  boundary  condition  will  be 
illustrated  for  the  condition 


o  (x,  0)  =  0  ,  x<  0  ; 

£*  Z 


(C3) 


the  procedure  is  analogous  for  the  remaining  boundary  conditions. 

In  light  of  the  foregoing  comments  regarding  interchange  of  differen¬ 
tiation  and  integration  in  (3.  32),  the  formula  for  o  ( x ,  z)  in  the 

zz 

interior  of  the  half-space  may  be  derived  from  (1.  4),  (3.  3),  (3.  32) 
and  (3.31)  as 
oo 

o  (x,  z)  =  f  E(M.)2ijtn  (K)(2H^-k^)(e  ns^K^Z-e  ndfK^z)e^KxdK, -o^<x<oo,  z>  0 

ZZ  S  S 

(C4) 

MA  k  (kf-k^)D+(k  ) 

E(h)  = -  0  R_  3  a - £ - —  .  (C5) 

2tt  ns(kR)ns(K)(K-kR)R(x)D"(K) 


oo 


where 


It  is  clear  that  °Z7,(X>  z )  as  given  in  (C4)  vanishes  at  z  =  0,  but  it 
must  be  proved  that  a__(x,z)  is  continuous  up  to  z  =  0.  Fix 

Z  Z  ni"*  1  — 

x  =xq<0  and  close  the  contour  in  the  lower  half  ?t- plane,  giving 
o  (x  ,  z)  =  [E(n)2ian  (x)(2ft2-k2)(e"ns^z-e  nd^H ^Z)e1KXOdx 

ZZ<  O  •'  s  s 

B- 

+  2iTi (Residue  at  -k^],  x<0,  z>0,  (C6) 

XV  o 

where  B_  is  the  branch-line  contour  in  the  lower  half  x-plane 
(see  Fig.  4).  The  residue  term  will  be  continuous  by  inspection  so  it 
remains  to  show  that  the  integral  in  (C6)  is  continuous.  The  contri¬ 
butions  to  this  integral  from  all  segments  of  B_  except  the  ones 
along  the  negative  imaginary  axis  represent  continuous  functions  of 
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z  since  they  are  integrals  aJong  finite  segments  of  a  function  contin¬ 
uous  in  both  z  and  k.  To  prove  continuity  of  the  contributions  from 
the  inmiite  segments,  first  change  the  variable  to  h  =  -in,  leaving 
two  integrals  in  which  n  runs  from  0  to  oo.  The  integrand  in  each 
will  be  0(rpe  o)  uniformly  in  z  as  r]-oo,  so  the  integrals  can 
easily  be  shown  by  the  M-test  to  be  uniformly  convergent  with  respect 
to  z  for  0<z<z^,  for  any  z^>0.  Now  apply  the  following  theorem 
(see  pp.  25-26  of  [24]): 

If  f(r),  z)  is  continuous  for  Osti^K,  0£z^zq  for  all  K>0,  and 
the  integral 

oo 

F(z)  =  Jf(nfz)dn 

o 

converges  uniformly  with  respect  to  z  for  0<z<zq,  then  F(z)  is  a 
continuous  function  of  z  in  this  interval,  the  continuity  being  one¬ 
sided  at  the  end-points. 

By  letting  f(r) ,  z)  equal  either  of  the  two  integrands,  F(z) 

then  represents  the  contribution  to  a  (x  ,  z)  from  the  respective 

zz  o  r 

segment  of  B_.  The  conditions  of  the  theorem  are  met  so  these 
contributions  are  continuous  up  to  z  =  0.  Thus  c  (x,  z)  is  continuous 
up  to  z  =  0  for  x<0,  the  limit  as  z-*0  may  be  taken  inside  the  inte¬ 
gral  in  (C4),  and  the  boundary  condition  (C3)  is  satisfied. 

The  functions  G+(?i)  and  F"(h)  are  defined  by  (3.  26)  and  their 
asymptotic  forms  are  obtained  in  (3.  27),  provided  that  order  condi¬ 
tions  (3.  2)  hold  as  x-*0  along  the  surface.  It  will  now  be  shown 
that  conditions  (3.  2)  hold  for  the  solution  given  by  (3.  3)  and  (3.  32). 
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When  (3.3)  and  (3.3  2)  are  substituted  into  (1.4)  and  the  inte¬ 
gration  contour  closed  by  large  arcs  in  the  upper  half  plane,  there 
results 


Oxz(x,  0)  =  2iri[X  Residues]  +  f  C(Jt)e*Hxdx, , 


B, 


9a 


(C7) 


xz 


9x 


z  =0 


=  ZiriC^  Residues]  +  Jix  C(K)eiKXdn  , 


B+ 


where  G(h)~P|x[^  as  |h|-oo  along  B_|_ ,  P  is  a  constant,  and  the 
integrands  are  continuous  along  B+.  The  residue  terms  and  the  con¬ 
tributions  to  the  integrals  from  all  segments  of  B^.  except  those 
along  the  positive  imaginary  axis  represent  functit  continuous  up 
to  x  =  0.  If  the  change  of  variable  h  =  ir)  is  made  in  the  portions  of 
the  branch  line  integrals  along  the  imaginary  axis,  they  take  the 

forms  oo 

Jc(H)e1KxdH  =  Jc^nJe^d-n  +  Od)  as  x-0,  (C8a) 


oo 

JiKC(K)e  Kxd»  =  -JriC':(ri)e‘Tlxdti  +  0(l)  as  x-0,  (C8b) 

B+  0 

■>'  ju  3 

where  C  (t))~P  r\2  as  T|-oo.  The  integral  in  (C8a)  is  clearly 
bounded  as  x-0.  To  analyze  the  integral  in  (C8b),  the  following 
extended  Abel  limit  theorem  may  be  used  (see  pp.  288-289  of  [26]): 

If  y  >  -  1  and  f(t)/AtY-l  as  t -oo,  then 


oo 


1  +  v  rf(t)e-3tdt 


Aliv+1) 


1  as  s  -  0*  . 
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When  this  is  applied  to  (C8b)  and  the  result  substituted  in  (C7),  it 


follows  that  8o  /9x 
xz 


_i  + 

-  G{x  2)  as  x-*0  . 


z~0 


An  analogous  procedure  may  be  used  to  verify  that 
=  0(1)  as  x-*  0  . 


a  (x,  0)+|iN9^u/8x^ 

XZ 


z=0 
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